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The gravitational wave signal from a binary neutron star inspiral contains information on the 
nuclear equation of state. This information is contained in a combination of the tidal polarizability 
parameters of the two neutron stars and is clearest in the late inspiral, just before merger. We use 
the recently defined tidal extension of the effective one-body formalism to construct a controlled 
analytical description of the frequency-domain phasing of neutron star inspirals up to merger. Ex- 
ploiting this analytical description we find that the tidal polarizability parameters of neutron stars 
can be measured by the advanced LIGO- Virgo detector network from gravitational wave signals 
having a reasonable signal-to-noise ratio of p = 16. This measurability result seems to hold for all 
the nuclear equations of state leading to a maximum mass larger than 1.97Mq. We also propose 
a promising new way of extracting information on the nuclear equation of state from a coherent 
analysis of an ensemble of gravitational wave observations of separate binary merger events. 

PACS numbers: 95.85.Sz, 04.30.Db, 04.40.Dg 



I. INTRODUCTION 

Binary neutron star (BNS) inspirals are among the 
most promising sources for the advanced version of the 
ground based gravitational wave (GW) detector network 
LIGO- Virgo. BNS's evolve under the influence of gravi- 
tational radiation reaction leading to a GW inspiral sig- 
nal whose amplitude increases up to the merger, while 
its frequency also increases up to a merger frequency 
/merger S 2000 Hz. One of the goals of the observation of 
GW signals from BNS systems is to improve our knowl- 
edge about neutron star (NS) structure and the highly 
uncertain equation of state (EOS) of NS matter. Ad- 
vanced LIGO is expected to be able to detect about 40 
BNS merger events per year [l] with signal to noise ratio 
(SNR) p > 8. The question that we shall address here 
is whether such observations can allow us to learn some- 
thing useful about the EOS of neutron star matter via 
the measurement of tidal polarizability parameters from 
the inspiral signal. 

In Newtonian gravity the (quadrupolar) tidal polariz- 
ability of a body is usually measured by means of the di- 
mensionless Love number &2 such that P2 = 2/(3G)fc2-R 5 , 
where R denotes the radius of the NS, yields the ratio be- 
tween the tidally induced quadrupole moment Q a b and 
the companion's perturbing tidal gradient G a b — d a bU. 
The generalization of the concept of tidal Love number 
&2 to strongly self-gravitating objects (NS or black holes) 
was discussed long ago by one of us as part of the the- 
ory of motion of compact bodies [2] . This work indicated 
how, by matching a quadrupolar deformed NS geome- 
try treated a la Thorne and Campolattaro [3] , one could 
compute ki for a given neutron star EOS. Recently, an 



explicit, simple, way of doing this matching computa- 
tion of &2 has been obtained by Hinderer [4]. The re- 
sulting numerical values for hi obtained in Ref. [4] have 
then been used in a preliminary analysis of the measur- 
ability of tidal effects in BNS GW inspiral signals [5]. 
However, this early work has been marred by a calcula- 
tional error in [4] leading to a substantial overestimate 
of the value of fe. Later work [6, 7] emphasized that 
&2 is a strongly decreasing function of the NS compact- 
ness C = GM/(c 2 R), such that k2(C) formally vanishes 
in the black hole limit 1 C — > 1/2, and generalized the 
computation of Love numbers so as to include gravito- 
magnetic tidal polarizability coefficients as well as higher 
multipolar contributions. Recently [9] the tidal polariza- 
tion parameter 2 p2 was computed for a wide range of 
EOS. Moreover, the question of discriminating between 
NS EOSs via GW observations with the advanced LIGO- 
Virgo detector network, using the early part (frequencies 
/ < 450 Hz) of the inspiral signal, has been discussed 
and answered in the negative in [9]: only if one has a 
GW signal with very high SNR p = 35, and if the actual 
EOS of NS matter is unusually stiff, can one start distin- 
guishing (at the 68% confidence level) the early-inspiral 
tidal signal from the noise. 



It was already mentioned in Ref. [2] that the &2 of a (four- 
dimensional) black hole vanishes. See [8] for a discussion of 
black- hole Love numbers in higher spacetime dimensions. 
We follow the notation for tidal polarizability parameters intro- 
duced some time ago in the General Relativistic Celestial Me- 
chanics formalism [10]: namely pn for the £ th -multipolar mass- 
type (gravito-electric) coefficient, and ag for the corresponding 
spin-type (gravito-magnctic) one. 



The reason why Refs. [5, 9] performed a conservative 
data analysis based only on the early inspiral GW sig- 
nal, / < 450 Hz, was that their analysis was based on 
using a purely post-Newtonian (PN) expanded descrip- 
tion of the phasing, without having any way of control- 
ling the validity of this description for frequencies above 
450 Hz. More precisely, they use a TaylorF2-type [11] 
description of the frequency-domain GW phase of the 
form #■(/) = #pn(/) +^n (/)> wrtn a point- mass phasing 
^pn(/) treated at 3.5PN accuracy [11] and with a tidal 
phasing &£(/) treated at leading, Newtonian order [5]. 

Recently, a new, improved description of the dynam- 
ics, waveform and phasing of compact binary systems has 
been developed based on the effective one body (EOB) 
formalism [12-17]. In particular, the way to extend the 
EOB formalism so as to include tidal effects has been 
presented in [18]. Let us recall that the EOB formalism 
is an analytical framework which combines several differ- 
ent theoretical results and approaches, and, in particular, 
contains resummed versions of the usual PN-expanded re- 
sults. Such a framework has proven to be a powerful tool 
for constructing analytic waveforms that agree with nu- 
merical simulations. In the binary black hole case, EOB 
waveforms arc in agreement with high-accuracy numeri- 
cal waveforms at the remarkable level of 0.01 rad up to 
merger [19-21]. In addition, the tidal-EOB formalism 
of [18] has been successfully compared to state-of-the- 
art numerical simulations of BNS systems [22, 23]. This 
comparison showed that the tidal-EOB formalism could 
reproduce the numerical phasing essentially up to merger 
within numerical uncertainties. 

This successful comparison (together with recent ana- 
lytical progress [24] in the computation of the EOB tidal 
interaction potential) motivates us here to use the tidal- 
EOB formalism as a way to define a controlled analyti- 
cal description of the phasing of tidally interacting BNS 
systems up to merger. More precisely, we will show be- 
low that the tidal contribution ^eob(/) *° the Fourier 
domain phase predicted by the tidal-EOB approach can 
be represented (within less than 0.3 rad) by a certain 
(PN-type) analytical expression up to merger. This will 
allow us to perform a data analysis using the full tidal 
phasing signal up to merger, while keeping the conve- 
nience of having an explicit analytical representation of 
the tidal phasing (instead of the well-defined, but more 
indirect, full EOB description of tidally interacting BNS 
systems) . Using such a EOB-controlled description of the 
tidal phasing up to merger, we will show (see Fig. 4 be- 
low) that the EOS-dependent tidal polarizability param- 
eters Gfj,2 of NSs can be measured, at the 95% confidence 
level, with the advanced LIGO- Virgo detector network 
using GW signals with reasonable SNRs (p = 16) for all 
EOS in the sample we shall consider (only restricted by 
the observational constraint of yielding a maximum mass 
larger than 1.97M Q [25]). In addition we shall propose a 
new way of extracting EOS-dependent information from 
a coherent analysis of a collection of GW observations 
of separate BNS merger events, which promises a large 



increase in measurement accuracy. 

In this paper we will focus on BNS systems, but the 
formalim we present can be used as it is for discussing the 
measurability of tidal parameters in mixed BH-NS binary 
systems. This would allow one to go beyond the recent 
works [26, 27] dealing with some aspects of the measura- 
bility of tidal polarizability coefficients from mixed binary 
systems. 

The paper is organized as follows: in Sec. II we will re- 
view the main elements of the tidal-EOB formalism and 
present the analytic tidal phasing model in the frequency 
domain that we will use in estimating the measurabil- 
ity of fi2- The theoretical aspects of our measurability 
analysis are given in Sec. III. The numerical results for 
the measurability of //2 are presented in Sec. IV, while 
concluding remarks are gathered in Sec. V. The paper 
is completed by two Appendices. In Appendix A we ex- 
tend and complete the review of the tidal-EOB formalism 
of Sec. II, giving in particular the explicit analytical ex- 
pressions for the tidal corrections to the EOB waveform. 
Finally, Appendix B collects the PN-expandcd formulas 
for the tidal phasing for a general rclativistic binary that 
are used in the main text. When convenient, we use ge- 
ometrized units with G = c = 1. 



II. ANALYTICAL TIDAL PHASING MODELS 
IN THE FREQUENCY DOMAIN 

The main aim of the present paper will be to estimate 
the measurability of tidal parameters by making use of 
the full BNS inspiral signal, including the late-inspiral 
part just before merger, where tidal effects are strongest. 
We will do so by taking advantage of the recent devel- 
opment of an analytical model which can accurately de- 
scribe the full inspiral signal up to merger. Indeed, in 
Refs. [22, 23] state-of-the-art numerical simulations of 
inspiralling BNS systems were compared to several an- 
alytical models. It was found that the EOB model (in 
its tidally extended version as defined in [18]) was able 
to match the numerical results up to merger. The EOB 
model (dynamics and waveform) is originally defined in 
the time domain. For the data-analysis purpose of the 
present paper it will be convenient to have in hands an 
analytic representation of the waveform in the frequency 
domain. The derivation of such an analytic frequency- 
domain phasing model will be the topic of the present 
Section. 



A. Tidal effects in EOB dynamics 

Let us recall that the EOB formalism [12, 13, 15] con- 
sists of three main elements: (i) a resummed Hamiltonian 
describing the conservative dynamics; (ii) a radiation- 
reaction force computed from the instantaneous angular 
momentum loss; (iii) a resummed waveform. 



For a nonspinning binary black hole (BBH) system of 
masses Ma, Mg the EOB Hamiltonian is given by 



H, 



EOB(r,p rt , Pv ) = Mc z Jl + 2v(H cB -l) , (1) 



where 



i + p j+*S 



(2) 



Here M = Ma + Mb is the total mass, v = 
Ma Mb /{Ma + Mb) 2 is the symmetric mass ratio, and 
z 3 = 2^(4 — Zv). In addition, we are using rescaled di- 
mensionless (effective) variables, namely r = r ab<? I GM 
and p v = P v c/(GMaMb), and p rt is canonically con- 
jugated to a "tortoise" modification of r [17]. The 
crucial input entering this Hamiltonian is the "radial 
potential" A(r), whose leading-order approximation is 
A(r) « 1 - 2/r + • • • = 1 - 2GM/(c 2 r AB ) + ■■■ 

The proposal of Ref. [18] for including dynamical tidal 
effects in the conservative part of the dynamics consists 
in using a tidally-augmented radial potential of the form 



A(r) =A u (r)+A tldal (r). 



(3) 



where A (r) is the point- mass potential defined in 
Eq. (A4) of Appendix A, while ^4 tldal ( r ) is a supplemen- 
tary "tidal contribution" describing the tidal interaction 
potential. In terms of the dimensionless gravitational po- 
tential u = GM/(c 2 rAB) = 1/r it reads 



\ tidal 



= £-' 

e>2 



Kju 2e+2 Af d *\ U ). 



(4) 



Here the term kJu 2+2 represents the multipolar tidal 
interaction of degree £, taken at Newtonian order in a 
PN expansion. The dimensionless EOB tidal parameter 
kJ entering Eq. (4) is related to the tidal polarizability 



coefficients G[ig ' of each neutron star as [18] 



where 



Kg 



{2£-l)\l 



Xb Gj-ig 

Xa~Jgm/^) 



2\2l+l ' 



(5) 



(6) 



where we recall that M — Ma + Mb denotes the total 
mass of the binary and Xa = Ma/M. The tidal polariz- 
ability coefficient Gfif has the dimension [length] 2f+1 . It 
measures the ratio between the £-th multipole moment 
induced in body A and the external tidal gradient felt 
by body A. Among these multipolar tidal polarizabil- 
ity coefficients, the dominant one is the quadrupolar, 
£ = 2, one, G[i£. [Note that [ii is denoted by A in 
Refs. [5, 9, 26]]. In addition, if Ra denotes the radius 
of body A, Gpif is related to the corresponding dimen- 
sionless Love number kg by 



(21-1)11 G[if 



2kg R A 



21+1 



(7) 



so that 



Kg 



2k, 



_ X_b ( Ra 

X A \G{M A + M B )/c' 



2l+\ 



(8) 



The additional factor A^ ldal (u) in Eq. (4) represents 
the effect of distance-dependent, higher-order relativistic 
contributions to the dynamical tidal interactions: 1PN, 
i.e. first order in u, 2PN, i.e. of order u 2 , etc. Here we 
will use the following "Taylor-expanded" form of A^ ldal 



\ tidal 



(it) = l + a^u + a^u 2 



(9) 



; .< f > 



where an' are functions of Ma, Ca, and kg for a general 
binary and are defined as (see Eq. (37) of [18]) 



-M = <H_ 



+ Kfa r , 



Kn T" K,i 



(10) 



AW 



where a n is the coefficient of the the nPN fractional 
correction to the tidal interaction potential of body A. 
(see Sec. IIIC of [18]). The individual dimensionless co- 
efficient a n is a function of the dimensionless ratio 
X A = M A /M. [Note that X B = M B /M = 1 - X A \- 
The analytical expression of the first post-Newtonian, 
quadrupolar (I = 2) coefficient a ± has been reported 
in [18] (and then confirmed in [28]) and reads 



A(2) 



= -X A - 
2 



(11) 



Recently, Ref. [24] has succeeded in computing the first 
post-Newtonian octupolar (£ = 3) coefficient a x , as 
well as the second post-Newtonian quadrupolar (£ = 2) 

and octupolar (£ — 3) coefficients a- . The most rele- 
vant 2PN quadrupolar coefficient reads 



AW 
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1 



= *r^ + 8^ + 3 - 



(12) 



In the equal-mass case, Xa = 1/2, the values of these 



coefficients arc a 1 ■ ' 

(2) 



O 



(2) 



5/4 = 1.25 and a 



A(2) 



a\ ' = 85/14 w 6.071429. A recent comparison [22, 23] 
between EOB predictions and BNS numerical simula- 

— (2) 

tions concluded that a 2 < 40. In the following, we shall 
restrict ourselves to considering only tidal quadrupolar 
contributions, i.e. we will take only the £ = 2 value 
in Eqs. (4) and (9). It is shown in Section A 2 of Ap- 
pendix A that the effect of higher-^ tidal corrections is 
small. It will be neglected in our analysis. 



B. EOB waveform and its stationary phase 
approximation 

When considering tidally interacting binary systems, 
one needs to augment the point-mass waveform h® m by 



tidal contributions. Similarly to the additive tidal mod- 
ification (4) of the A potential, we will here consider an 
additive modification of the waveform, having the struc- 
ture 



him — he 



' tidal 



(13) 



See Appendix A for the explicit expressions of h® m and 
hfa ■ In turn, this tidally modified waveform defines 
a corresponding tidally modified radiation reaction force 
T v through its instantaneous angular momentum loss. 

The radiation-driven EOB dynamics defined by 
-Heob(^4) and T v (where both A and F v are tidally 
modified) allows us to compute a time-domain multi- 
polar GW signal h lm (t) = A lm (t)e-' l ' t > t ™ { - t \ Following 
Refs. [22, 23], we characterize the (time-domain) phas- 
ing of the quadrupolar waveform hi2 {t) by means of the 
following function of the instantaneous quadrupolar GW 
frequency uj = uj(t) = d<p/dt [where </>(£) = 022 (t)] 



M 



dm _ {low 



dliiuj(t) uj(t) 



(14) 



In the stationary phase approximation (SPA), the phase 
&(f) of the frequency-domain waveform /i(/), i.e. the 
phase of the Fourier transform of the time-domain 
(quadrupolar) waveform, 



M/) = A(f)e-^f\ 



(15) 



is simply the Lcgendrc transform of the quadrupolar 
time-domain phase 0(t), namely 



#spa(/) = 2nft f - 4>it f ) 



IT 

4' 



(16) 



where tf is the saddle point of the Fourier transform, i.e. 
the solution of the equation w(£/) = 27r/. Differentiating 
Eq. (16) twice with respect to / leads to the following 
link between ^spa(/) and the function Q u (oS) 



d 2 V SPA (LU f ) _ Q u (u f ) 



dhj 1 . 



(17) 



7 



where uif now denotes the Fourier domain circular fre- 
quency uif = 2?tf. Below we will simply denote the 
Fourier domain frequency w/ as lo without bothering to 
distinguish it from the time-domain co(t). 

In the following we shall decompose the result (17) in 
its point-mass and tidal parts, thereby relating the "tidal 
part" of the Fourier-domain phase ^spa(/) to the "tidal 
part" of <9cj(w). On the one hand, the tidal part, say 
Q^(w), of Qu,(ui) is computed as 



Q«(w) = Q«(w) - 0°(w), 



(18) 



where Q^(uj) is the outcome of a point-mass EOB simu- 
lation, i.e., one without tidal effects in both the dynamics 
and the waveform. Then, the corresponding tidal part, 

*bob(«) = ^eobH - S'eobH, (19) 



of the Fourier-domain EOB phase satisfies, within the 
SPA approximation, the relation 



rf2 ^EOB SPA 



H CM 



dcu 2 



(20) 



Let us emphasize that we expect the SPA approximation 
to the phasing to remain accurate up to the merger. In- 
deed the small parameter that controls the validity of the 
SPA is essentially e a di a b = w/oj 2 = 1/Q U . For instance, 
Ref. [29], Eqs. (3.9)-(3.10), has computed the next con- 
tribution beyond the leading SPA and found that it in- 
troduces a dephasing 5W which, in the case of Newtonian 
chirps, is equal to 5V = (23/24) x (2/9) x 1/Q U . The 
quantity Qa, is very large during early inspiral and de- 
creases towards the merger. Looking at the value of the 
full Qu = Q° + Q^ m the exact EOB description of 
tidally interacting BNS systems, we have checked that 



the (equal-mass) value of Quj(uj) for uj 



(where 



^contact j g ^ ne eob approximation to the merger fre- 
quency, see below) remains larger than about 20 for real- 
istic compactnesses (C = 0.14 — 0.18). Though this value 
is reduced (by ~ 10) from the corresponding point-mass 
value Qji(a; con tact), it is still comfortably large compared 
to 1, so that one can expect the phasing error linked to 
the use of the SPA to be a small fraction of a radian. 



C. PN-expansion of the EOB tidal phasing 

In Sec. IV below we shall estimate the measurability of 
the tidal parameter k^ by computing the Fisher matrix 
F corresponding to the simultaneous measurement of a 
tidal parameter, say At oc k% , with several other, non 
tidal, parameters, say A a , a = 1, . . . , n. Though, in prin- 
ciple we could numerically compute the relevant Fisher 
matrix F by evaluating the numerical derivatives of the 
full, Fourier domain, EOB waveform /ieob(/; At, A a ) 
with respect to all the parameters (At, A q ), it will be con- 
venient to estimate F by replacing the numerically com- 
puted /ieob(/; AT,A a ) (which involves computing the 
numerical Fourier transform of a numerically generated 
time-domain EOB waveform /ieob(£; At,A q )) by some 
sufficiently accurate analytic approximation. We will do 
so by combining several approximations, the validity of 
which we shall control. The first approximation we shall 
use is the SPA, which we have discussed in the previous 
section. The second approximation will consist in using 
post-Newtonian expansions to derive adequately accu- 
rate expressions of the two parts of the Fourier domain 
phase 



V(f)=*°(f:\ a )+* T (f; A T ,A„). 



(21) 



In this section we study how many terms in the PN 
expansion of the tidal phase W (/; At, A ) we must re- 
tain to obtain an approximation to tf' T (/) which remains 
reasonably close to the EOB prediction up to merger. 
From Eq. (19) we see (in the SPA) that to answer this 



question we need to compare the PN-expansion of the 
tidal part of Q^(w) of Q u (u>) to the "exact" value of 
Qu(w) defined by the EOB model. During most of the 
inspiral , not only is the phase evolution quasi-adiabatic, 
i.e. Q u ^> 1, as already discussed above, but the dy- 
namical evolution can also be well approximated by an 
adiabatic quasi-circular inspiral. In the latter approxi- 
mation, the function Q^ = ur /Co is obtained by writ- 
ing the balance equation between the instantaneous en- 
ergy flux F at infinity and the adiabatic evolution of the 
energy of the system (i.e., the Hamiltonian, H(ui)) ex- 
pressed as a function of the instantaneous GW frequency 
u = 2 J? (where fi is the orbital frequency). This yields 
— F = dH/dt = (dH/dco)Co, from which one obtains 



-jadiab/, ,\ _ , .2 



QT a » 



dH{co)/dio 



(22) 



Reexpressing this result in terms of the dimensionlcss 
rcscalcd angular momentum j = J/(G/j,M), the Newton 
normalized energy flux F = F/F Nowton = F v /T™ cwton , 
and replacing the independent variable u by the usual, 
dimensionless PN ordering parameter 



lGMu\ 



2/3 



ttGM/ 



2/3 



(23) 



leads to an expression of the form 



Q adiab (x; v, A T ) = — x" 5 / 2 6(x; v, A T ) (24) 

40^ 



where the function b(x\ v, At), defined as 



b(x; v, Xq 



-2x 



F(x; v, At) 



(25) 



is simply equal to 1 in the Newtonian approximation. 
More precisely, it starts as 

b(x; u, At) = 1 + — (743 + 924^) x - 47ra; 3/2 + 0{x 2 ). 
336 

(26) 

Starting from the adiabatic EOB dynamics, the func- 
tion j(x) is obtained by eliminating u = \/r between 
the EOB expression j 2 (u) = — A' (u) / (u 2 A(u))' (obtained 
by minimizing the effective potential for circular orbits 
A(r)(l + j 2 /r 2 )) and the expression of Q in terms of u 
obtained from the Hamilton equation ft = dH/dJ .[See 
Sec. Ill and IV of Rcf. [18] and Appendix A for more 
details about the EOB circular dynamics] . On the other 
hand, the function F(x) is obtained by as a sum of var- 
ious resummed circular multipolar waveforms h£ m (x) of 
Rcf. [16]. 

In the following, we shall replace the (EOB-rcsummed) 
adiabatic approximation Q^ diabatic , Eq. (22), to Q^ OB , by 
a sufficiently accurate PN expansion of Q adlabatlc . This 
is to avoid an inaccurate feature of Q adlabatlc during the 
late inspiral. From Eq. (22) we see that Q adiabatic ( w ) is 



proportional to the derivative dH(uj)/duj which, by con- 
struction, vanishes at the Last Stable Orbit (LSO), where 
the circular energy H(ui) reaches a minimum. By con- 
trast, the exact Q^ OB does not vanish at the LSO, nor 
the PN expanded version of Q^ OB that we shall use. The 
frequency corresponding to the tidal-EOB defined LSO 
happens to be quite close to the contact frequency. Using 
Qadiabatic U p -j- Q corL t ac t might then introduce inaccura- 
cies in the phasing just before merger. Our use below of 
a suitable PN-expanded representation of Q M avoids this 
source of uncertainty and maintains consistency with the 
SPA by allowing the value of Q^ at contact to remain of 
order 20 for all cases considered. 

Current analytical knowledge that has been incorpo- 
rated in the EOB description of tidal effects [22, 23] al- 
lows us to compute the tidal part Q^{uj) of Q u (ui) and 
therefore, using Eq. (20), the tidal part \P T {uS) of the 
Fourier domain phase ^(w) beyond the 1PN accuracy 
obtained in Ref. [30]. First, the fact that the EOB for- 
malism naturally accomodates the inclusion of tail effects 
in the waveform allows us to obtain a PN-expanded tidal 
phasing model that is analytically complete up to 1.5PN 
order. [In addition, the EOB formalism already contains 
the next order tail effects at 2.5PN order]. Second, the 
EOB approach is designed in a way which makes it easy 
to complete it beyond current analytical knowledge by us- 
ing effective field theory methods. In particular, Ref. [24] 
recently computed the 2PN tidal contributions to the 
EOB radial potential A(u), i.e. the coefficient «2 of u 2 
in Eq. (9) (see Eq. (A10)). As mentioned in Ref. [24], 
a straightforward extension of the method used to de- 
rive the 2PN tidal contribution to A(u) can allow one to 
derive the 2PN tidal contribution to the waveform. How- 
ever this calculation has not yet been completed. Waiting 
for this result, we shall here use the natural flexibility of 
the EOB formalism, to parametrize the 2PN tidal cor- 
rections to the multipolar waveform by means of some 
parameters that we will call /3^ m . Let us recall that in 
order to obtain Q^ at, say, the fractional 2.5PN accuracy, 
the energy flux F must be computed by retaining all the 
£ = 2 and £ = 3 multipolar contributions to the wave- 
form. Then, to obtain the flux to 2.5PN accuracy we need 
the quadrupolar £ = in = 2 waveform (stripped of its tail 
factor) to 2PN fractional accuracy, and the odd-parity 
i = 2, m = 1 and £ = m = 3 and £ = 3, m = 1 even- 
parity waveforms at 1PN fractional accuracy. Following 
Refs. [18, 23], we shall parametrize such highcr-PN tidal 
corrections to the waveform along the following model 



tidal 



with 



h 



Atidalfjowt 



, A tidal 



B tidal 



(27) 



hf^ 1 = ^ 1MI »~' (1 + p( m (X A )x + P l 2 m {X A )x 2 + ...) 

(28) 
For the time being, the only such PN fractional tidal 
waveform correction which is known is P 22 (Xa)- Using 
the lPN-accuratc results of Rcf. [30] one can indeed dc- 



rive the following explicit analytical expression 

22 -202 + 560X A -340A 2 +45Al 
A ™ = 42(3-2^) ' (29) 

However, at 2.5PN, the final result depends on sev- 
eral other higher-corrections, namely (3 2 2 (X A ), /? 21 (^Gt), 
Pi 3 (Xa) and /3 31 (X A ) , that account respectively for 2PN 
fractional tidal corrections to the t = m = 2 multipolc 
and for 1PN fractional tidal corrections to the I = 2, 
to = 1 and £ = 3, m = 1 and m = 3 subdominant mul- 
tipolcs. In Appendix B we present the explicit expres- 
sions for Q^ and & T at 2.5PN accuracy for the general 
case of unequal mass binary systems. In the text below 
we shall specify those general formulas to the particular, 
but physically most relevant, case of equal-mass neutron 
star binaries (having therefore equal-compactnesses and 
equal tidal parameters). 

In the equal-mass case, because of symmetry reasons, 
the only higher-order tidal waveform parameter that con- 
tributes to the phasing is /3 22 . We then arrive at the 
following explicit expression for the 2.5PN accurate tidal 
contribution to Q^ 
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where we recall that a y 22 > = a 2 K "'{X A = 1/2) = 85/14 



and where /3 22 denotes the value of the (unknown) func- 
tion Pl 2 {X A ) for X A = 1/2. Using Eq. (20), the corre- 
sponding 2.5PN accurate tidal phase of the Fourier trans- 
form of the GW signal reads (in the SPA) 
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Such an explicit representation of a Fourier domain phase 
as a polynomial in x(f) ex J 2 / 3 is usually called Tay- 
lorF2 [11]. The 2.5PN TaylorF2 formula (31) improves 
the 1PN result of Ref. [30] in that: (i) tail effects are 
included up to 2.5PN order; (ii) a large part of the 2PN 
term is explicitly computed, although it still depends on 
the yet uncalculated quantity j3 2 2 (2PN tidal correction 
to the waveform). Note that at leading, Newtonian, or- 
der, Eq. (31) predicts that the (equal-mass) tidal dephas- 
ing at contact, i.e. for ^contact = Ca — Cb (see Eq. (36)) 
is of order 
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which, for the typical values k 2 = 0.08 and C = 1/6, 
yields —8.6 rad. With the further amplification of PN ef- 
fects discussed below this means that the tidal dephasing 
at contact is of order —10 rad (see Fig. 1). 

Let us now indicate why we expect that the contri- 
bution to the tidal phase coming from /3 22 is likely to 
be numerically subdominant compared to the currently 
known terms. Let us first note that at leading, Newto- 
nian order the overall coefficient (39/4)^ in the tidal 
phase Eq. (31), is, in view of Eq. (22), the sum of a tidal 
contribution from the Hamiltonian H and a tidal contri- 
bution from the energy flux F. More precisely, one finds 
that 



39 



(9 ff + 4 F ) 



(33) 



where the indices indicate the origin (H or F) of the con- 
tribution. Already at this leading-order level, one notices 
that the contribution from the energy flux is subdomi- 
nant (by a factor 2.25) with respect to the contribution 
from the Hamiltonian, i.e. from the radial potential A(u). 
When pursuing this analysis at the 1PN level and con- 
sidering the fractional PN modification of the tidal phase 
^2 T 5PN = ^2 T 5PN/^Newt = 1 + 0(x) in Eq. (31) one finds 
(still for the equal-mass case) 
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where we decomposed the 1PN fractional contribution 
into three parts: i) one coming from the leading order 
tidal terms in H and F; ii) one coming from the 1PN tidal 
correction to H (term oc a\ ); and (iii) one coming from 
the 1PN tidal correction to the quadrupolar waveform 
(and flux, term ex /3 22 ). In the equal-mass case one has 
a[ 2) = 5/4 = 1.25 and /3 22 = -11/672 w -0.0164. As a 
consequence of these numerical values, we see that: a) the 
coefficient of /3 22 is smaller by a factor 3.2 than the coeffi- 
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cient of a\ ; b) in addition, as the numerical value of /3 2 
is s=s —0.0164, its contribution to the total 1PN fractional 
coefficient is 0.10989 x 0.0164/2.05709 w 8.8 x 10~ 4 times 
smaller than the first term 2.057 and 7.2 x 10~ 4 times 
smaller than the sum of the first two contributions. Per- 
forming a similar analysis at the 2.5PN level (now insert- 
ing the known numerical values of a\ and /3 22 ) yields 
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Again we see that the contribution from /3 22 is likely to 
be subdominant. Indeed, not only is the coefficient of 
/3| 2 4.3 times smaller than the one of a 2 , but it is also 
about 149 times smaller than the known 2PN coefficient 



6.99 + 0.255 



(2) 



.51. Independently of these numerical 
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FIG. 1. Successive PN approximants to the (SPA) tidal con- 
tribution to the EOB phase as function of the £ — m = 2 
GW frequency. The leftmost vertical line indicates 450 Hz 
for a 1.4M + 1.4M© BNS system. The rightmost vertical 
line indicates the contact frequency that is taken as a fiducial 
analytical definition of the moment of merger. The plot refers 
to a C — 0.16, 7 = 2 polytropic model. 



ity. The phase ^eOB S pa( w ) * s com P u ted by integrating 
numerically Eq. (20) starting from the frequency ujq that 
marks the beginning of the inspiral waveform obtained 
when solving the EOB equations of motion numerically. 
This integration is done using the 2.5PN result for & T 
and <N/ T /du> as initial boundary conditions, and thus ojq 
needs to be chosen sufficiently small (i.e., the EOB in- 
spiral waveform has to be sufficiently long) so to have 
Q2.5PN w qEOB_ We re f er the reader to Appendix A 1 to 
get further technical details related to Fig. 1. 

The various PN approximations gathered in Fig. 1 are 
obtained from Eq. (31) and are represented as: thick 
dashed line (black online) for the Newtonian, dotted line 
(blue online) for the 1PN, dash-dotted line (red online) 
for the 1.5PN, dashed line (red online) for the 2PN and 
solid line (red online) for the full, 2.5PN phase. The 
leftmost vertical line indicates the frequency 450 Hz (used 
as cutoff in Rcfs. [5, 9]), while the rightmost vertical line 
indicates the frequency of "bare" contact, that defines 
within the EOB formalism the merger frequency. This 
bare contact is defined as the GW frequency where the 
relative distance R = M/u ps M/x is equal to the sum 
of the radii of the two NS R a + Rb = M a /C a + M B /C B , 



arguments, let us note that, as already mentioned above, 
an important feature of the adiabatic approximation to 
the function Q u (lu) is that it vanishes at the adiabatic 
LSO. Since tidal effects strongly influence the LSO loca- 
tion (see Ref. [18]) this indicates that tidal corrections 
to the Hamiltonian (i.e., to A(u)) have a dominant influ- 
ence on the shape of the Qu(uj) function below the LSO 
frequency and thereby on the tidal corrections to !^ T (w). 
In view of these arguments, in the following we will ne- 
glect the effect of /3| 2 both in the exact EOB phase and 
its 2.5PN approximant, ^spn^)- We will then work 
with Eq. (31) with j3f = (but af ] = 85/14 rs 6) as a 
numerically acceptable approximation to the 2.5PN tidal 
phase. 



D. Accuracy of PN-expanded representations of 
the EOB phasing 

Let us now study to what extent the "exact" EOB 
tidal phasing ^eOB spa obtained by integrating Eq. (20) 
with the exact (time-domain) Qu,(uj) defined by Eq. (14) 
on the right hand side, can be represented by various PN 
expansions. Fig. 1 illustrates the performance of Eq. (31) 
(considered at different orders of truncation) in reproduc- 



ing #j 



EOB s 



The figure refers to an equal-mass binary, 
with 7 = 2 polytropic EOS p = K/4 7 , with compactness 
C = 0.16. In addition, each star has M A = 1.47838Af Q , 
R A = 13.6438 km, k£ = 0.063122 so that the EOB di- 
mcnsionless tidal parameter of the system is equal to 
kT> = 75.2476, which yields Gjjl 2 = 19896km 5 . As men- 
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from which the gravitational wave frequency at con- 
tact is computed using Mw con t a ct = 27rAf/ contact = 
2(a; C ontact) 3 ^ 2 - In the equal-mass case Eq. (36) yields the 
simple result x contact =C A =C B . 

Among the useful informations contained in this fig- 
ure let us note that: i) the Newtonian approximation 
substantially differs from the EOB phase even at low fre- 
quencies, and exhibits a discrepancy of about 3 rad at 
merger; ii) as the PN order n is increased, the conver- 
gence towards the EOB prediction is non monotonic and 
the sign of the difference A n = !f n pN - $eob alternates as 
n takes the successive values 0, 1, 1.5, 2. This is linked to 
the alternating signs in ^spn m Eq. (31). In particular 
the difference A\ reaches the value —0.6 rad at contact; 
iii) it is only at 2.5PN accuracy that we get a rather ac- 
curate representation of the EOB tidal phase. Note that 
at merger, where the frequency parameter x reaches the 



value 



= C = 0.16, the fractional PN modification 



of the tidal phase is equal to 
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(37) 
which illustrates the effect of the successive PN approx- 
imations, labelled here by the corresponding PN order, 
(1,1.5,2,2.5). 

To firm up the conclusions drawn from the particular 
model (with compactness C — 0.16) considered in Fig. 1, 
we studied what happens when the compactness varies 
within a realistic range. Let us recall that the magnitude 




FIG. 2. Difference between the 2.5PN-expanded tidal 
(Fourier) phase and the corresponding exact EOB one ob- 
tained by integrating Eq. (20). Each curve refers to a 7 = 2 
polytropic model with different compactness. The vertical 
lines indicate the corresponding contact frequency. 



of the dimensionless EOB tidal parameter n^ = k^ + k§ 
is related (in the equal-mass case) to the Love number fe 
and to the compactness C by 



8C 5 



(38) 



[For a general multipolar index one has kJ = 
l/(2 2i ~ 1 )k e /C 2i+1 }. For a given EOS, as C increases, 
&2 decreases in a correlated manner [6, 9], so that 
K2 varies by about a factor 9 in a range of realis- 
tic compactnesses. For instance, in the case of the 
7 = 2 polytrope that we are currently discussing, as C 
varies between 0.14 and 0.19, n^ decreases from 183.37 
down to 21.757, with radii correspondingly varying from 
14.369 km down to 12.435 km. We generalized the com- 
parison reported in Fig. 1 for three different compact- 
nesses C = {0.14,0.16,0.18}. The results for the differ- 
ences 

Z\™ EOB <FH = ^ 5 pnH - 8£ OB ( w ) (39) 



III. MEASURABILITY OF TIDAL 
PARAMETERS: THEORETICAL DISCUSSION 

A. Fisher matrix formalism 

Under usual simplifying assumptions (Gaussian noise, 
sufficiently high SNR) the variance a 2 T in the mcasurc- 

ment of k^ is computed using the standard Fisher ma- 
trix formalism, as already used in the context of binary 
systems in Refs. [5, 9, 26, 31, 32]. When considering a 
Fourier-domain waveform h(f; Xi) that is a function of 
?i + l parameters {A.;} = {A a , At}, the Fisher information 
matrix is a (n + 1) x (n + 1) matrix whose elements are 
given by 

Fa = (d Xi h,d Xj h) (40) 

where (, ) denotes the Wiener scalar product between two 
signal h(f) and k(f), denned as 
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with S n (f) denoting the one-sided strain noise of the de- 
tector. In absence of specific prior, the variance in the 
measurement of each parameter Xi is given by the corre- 
sponding diagonal element of the inverse Fisher matrix 
(or covariance matrix), 
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Assuming the SPA approximation and neglecting rela- 
tivistic corrections to the amplitude, the Fourier trans- 
form of the waveform is 
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where 6*(/ contact — /) denotes a Heaviside step function in- 
dicating that we cut off the inspiral signal above the con- 
tact frequency defined above in Eq. (36). [ This is a coarse 
approximation to the post-merger signal that might be 
improved by extending the EOB representation to an ef- 
fective description of the post-contact GW signal] . The 
amplitude parameter A has been shown to be uncorre- 
cted to the other parameters [31, 32], so that we shall 
forget about it in the following 3 . 

From this equation, the squared signal to noise ratio 
(SNR) is written as 
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are shown in Fig. 2. In all cases these differences are 
rather small, being < 0.3 rad at merger. In addition, let 
us note that the positive sign of aP^eob^t (gj ven ^g 
fact that W T is negative) means that using &25PN instead 
of the more exact EOB phasing S'eob * s a conservative 
way of estimating the mcasurability of tidal parameters. 



3 However, strictly speaking, because of the dependence of the 
step-function g(j contact —fj ; n Eq. (43) on t ne dynamical param- 
eters via / contact (A a ), there will be a small correlation between 
the amplitude parameter A and the other parameters. Follow- 
ing [31], we shall neglect this correlation which is not expected 
to modify our conclusions in any significant way. 



where f c = jc 011 ' 3 - ^ and the elements of the Fisher ma- 
trix are 
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In view of the common proportionality of both p 2 and i 7 ^ 
to the (randomly vaying) squared signal amplitude A 2 , 
it is convenient to define the following "reduced" Fisher 
matrix 
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One then sees that this reduced Fisher matrix can be 
written as 



(47) 



Fij= /4f7(/)0A,#(/)0A i #(/). 



in terms of which we can write the rescaled Fisher matrix 

as 



Pij = (d Xi *\d x .&) 



(50) 



The elements of the inverse of this reduced Fisher matrix 
then give the "SNR-normalizcd probable errors" , p a Xi on 
each parameter, namely 



a Xi = pa x% = \I{F 



(51) 



In the following S n is taken to be the ZERO_DET_tiigh_P 
anticipated sensitivity curve of Advanced Ligo [33] . The 
minimum of the effective dimensionlcss strain noise 
fS n (f) for this sensitivity curve is located at frequency 
/o = 56.56 Hz. In the following we shall often work with 
the reduced frequency parameter / = ///o- 



B. Phasing model and parameter dependence 
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Note that this measure is normalized to unity, 

r 

J j(f)df = 1. This measure naturally leads to defin- 

o 

ing a new (Euclidean) scalar product among real phase 

functions 



r 

(a\b) = JdMf)a(f)b(f), 
o 



(49) 



Concretely, we shall use a Fourier domain waveform of 
the type of Eq. (43), with a phase &(f) in the form 

nf) = *\f)+* T {f), (52) 

where, as above, &°(f) denotes the point- mass contri- 
bution to the SPA phase and W T {f) the tidal part. We 
shall approximate both contributions with some PN ex- 
pansion. As already discussed above, the tidal contribu- 
tion will be approximated by the 2.5PN accurate expres- 
sion of Eq. (31). Concerning the point-mass phase, it is 
currently analytically known up to 3.5PN order [34, 35]. 
As we shall further discuss below, for the purpose of 
the present paper it will be enough to use the follow- 
ing 2PN [36] accurate representation of the point-mass 
phase 
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where 
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and where the parameters the (Aj)[i..4i have the following 
meaning 



Ai = -<f> c - 7r/4, 

A2 = fot c , 

X 3 =M, 

A4 = V. 



(55) 
(56) 

(57) 
(58) 



Here <p c is a reference phase and t c a reference time, 
(54) M = ^/5jf i s the chirp mass and v = M A M B /M 2 
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the symmetric mass ratio. In addition, the parameter ft 
is a spin-orbit parameter and a a spin-spin one. 

As for the (quadrupolar) tidal contribution & T (f) we 
can write it in various forms depending on the choice 
of tidal parameters we want to fit for. For instance, if 
we choose as tidal parameter At determining the overall 
scale of the tidal phase the following symmetric com- 
bination of the two £ = 2 tidal polarizability coeffi- 
cients G\i 2 ' 
[length] 5 ), 



(2/3)fc 2 ' R\ b (with the dimensions of 
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we obtain a tidal signal of the form 
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pressed as a function of the chirp mass A3 and the sym- 
metric mass ratio A4. In this form the A a -dcpcndcncc 
of the factor v^ 3 A4 /(A/f A4 A,i) introduces correlations 
when fitting for At together with the A a 's. An alter- 
native choice might be to consider as tidal parameter the 
(dimensionless) combination 



unequal-mass case discussed below) we shall neglect the 
variability of the A a 's entering the PN-correction factor 

^TspnC/i ^o)' (i- c - ne gl ec t their small contribution to 
d\ a \fr T computing the Fisher matrix). 

In this work we keep in the tidal signal ^J 5 pn(/) only 
the contribution associated to the quadrupolar tidal de- 
formation as measured by k 2 or Gfif ■ Actually, the EOB 
formalism takes into account higher multipolar tidal in- 
teractions, as already done in previous work [18, 22, 23]. 
Using this theoretical result, we show however in Ap- 
pendix A 2 that the numerical contribution of higher mul- 
tipolc moments (£ = 3, 4) to the tidal signal is rather 
small {A\P T ~ —0.2 rad), so that we are entitled to ne- 
glect it to estimate the measurability of k 2 . However, 
we recommend that in fitting real GW signals to tidal 
EOB templates one includes also the higher multipolar 
tidal contributions. But, in order not to introduce new 
v \ a /^^JsPNt/^a)' parameters to be fitted, one should express the higher- 
order polarizability parameters, Gfie, in terms of G\i 2 
only. More precisely, using G\i 2 ~ k 2 R 5 , G113 ~ k^R 7 
and Gfi4 ~ k^R? , one can reexpress G/X3 and G/^4 in 
terms of G\x 2 and of the following combination of Love 
numbers 
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so as to minimize the correlations when fitting A T to- 
gether with the A a 's. Note, however, that there will al- 
ways remain correlations due to the A a -dcpendcnce in 
the fractional PN correction factor 'F 2 T 5PN (f;X a ). The 

A a -dependence of ^JspnI/; ^a) has several sources: i) 
an explicit dependence on M and v coming through the 
argument v = «a 3 ,a 4 / ; ii) an implicit dependence on 
the mass ratio coming from the individual Love numbers 
and the radii entering the definition of ^spn ( see Ap- 
pendix B). However, when computing the corresponding 
Fisher-matrix elements involving the partial derivative 
with respect to A a of &, the contributions from the Aa- 
dependence of the tidal part <P T are largely subdominant 
compared to the large, early-inspiral dominating contri- 
bution coming from dx^ . We have indeed checked 
that taking into account the variability of the A a 's within 
^25Pn(/' ^o) or neglecting it only changes the error <j\ T 
at the 1.5 x 10~ 3 fractional level. As for the A a variability 
of the prefactor of ^J 5PN in Eq. (60), when using (59) as 
tidal parameter, it was found (when ft and a are fixed, or 
well constrained) , because of the signs of the correlations 
between At and A a , to lead to a small, 0(2%), improve- 
ment in the measurability of At compared to that of A T 
given by Eq. (61). In the following we shall fit for At, 
Eq. (59), taking into account the variability of the A a 's in 
the prefactor of the tidal phase [as displayed in Eq. (60)], 
however, for simplicity (and easier comparison with the 



e.g., G/i3 ~ fc 3 / 2 (G/Z2) 7 / 5 . Replacing then the modified 
Love numbers k^/ 2 and £4/2 by some constant numbers 
(say of order 0.7 so as to approximately mimick the re- 
sult of realistic EOSs) we end up with an approximate 
description of higher multipolar contributions that is en- 
tirely expressed in terms of the quadrupolc polarizability 
parameter G/i 2 . 

Let us now comment on the form, Eq. (53), of the 
point-mass phase that we shall use in this work. This 
point-mass phase is only 2PN accurate. The reason 
for limiting our accuracy to this level is that, as we 
will sec explicitly below, the terms in the Fisher ma- 
trix that determine the measurability of the two dynam- 
ical parameters entering the point-mass phase, namely 
the chirp mass A3 = M. and the symmetric mass ra- 
tio A4 = v, are essentially 4 proportional to integrals of 
the following types: 7_io = f d\nf fj(f)v(f)~ 10 and 
7_6 = fdlnffj(f)v(f)~ (i . While the integral giving 
the signal to noise ratio, Eq. (44) is proportional to 
la = J dlnffj(f) and is roughly concentrated around 

a couple of frequency octaves around / = f/fo = 1, the 
integrals i_io and I-§ are mainly concentrated towards 
(different) lower frequencies. The concentration on the 



4 Here, for illustrative purposes, we keep only the leading-order PN 
signal contributing to the corresponding Fisher matrix element: 



e.g. dx^ ~ V 5 leading to ^33 ~ (y 5 \v 5 ) = I- 
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FIG. 3. Integrands, per frequency octave, of the inte- 
grals determining the measurability of M, v, p (SNR) and 
At. While most of the SNR is gathered around frequencies 
/ = //(56.56 Hz) ~ 1, the measurability of M and v is con- 
centrated towards lower frequencies (/ = f/fo < 1), and that 
of the tidal parameter At gets its largest contribution from 
the late inspiral up to the merger. The rightmost vertical line 
indicates the merger frequency for C — 0.1645, while the left- 
most vertical line marks 450 Hz for a 1.4Af Q + 1.4A/ Q BNS 
system. 



logarithmic frequency axis of several relevant measura- 
bility signals is illustrated in Fig. 3. Note in particular 
how the integrands of l-io (chirp mass) and I-q (sym- 
metric mass ratio) arc peaked at frequencies below the 
SNR integrand of Iq . Physically, this corresponds to say- 
ing that most of the useful cycles for the measurability of 
M. and v come from the early inspiral. As the PN expan- 
sion converges reasonably well for such low frequencies, 
using a 2PN accurate phasing is guaranteed to be a rea- 
sonably good approximation for the point-mass part of 
the phase. This has been checked by Ref. [32] for the 
measurement of M. and v, which found (see their Ta- 
ble II) that using a 2PN accurate (instead of a 1.5PN 
accurate, as in Ref. [31]) template led to only ~ 10% 
differences in the fractional uncertainties in v and M. 
We found, as expected, that the situation is even better 
for the measurement of At: namely, we found that the 
fractional uncertainty on At is changed (and actually im- 
proved) when using a 2PN template for if -0 , rather than 
a 1.5PN one, only at the 5 x 10~ 3 level By contrast to 
the cases of M. and v, the measurability of the tidal pa- 
rameter At is associated in the Fisher matrix to an in- 
tegral of the type I +10 = / d\nff^(f)v(J) w , which gets 
its largest contribution from the late inspiral up to the 
merger (see solid line in Fig. 3). More specifically, the in- 
tegrand of J+io, i.e. ex /7(/)/ 10/3 is equal to f 2 /S„(f). 
The ZER0_DET_high_P advanced LIGO noise curve S„(f) 
happens to be a rather flat function of / between ~ 50 Hz 



and ~ 800 Hz and then increases to reach a shot noise 
behavior S n (f) oc f 2 at high frequencies. This implies 
that the integrand of i+io, i.e. f 2 /S n (f), roughly grows 
like f 2 between 50 Hz and 800 Hz, to then asymptote 
towards a finite limit at high frequencies. The clear sep- 
aration between, on the one hand, the two SNR curves 
associated to M. and v (which are relatively close to each 
other) and on the other hand the SNR curve associated to 
At also indicates (as we shall discuss below) that M and 
v are strongly correlated among themselves, while At is 
not so strongly correlated to M and v. The figure also 
displays two possible cut-off frequencies for the measure- 
ments of the tidal signal: the conservative value 450 Hz 
(dashed vertical line) used in Refs. [5, 9], / = 7.956, or 
the compactness-dependent contact frequency that we 
shall use here, (7rM/) contac t = C 3 / 2 (dash-dotted verti- 
cal line, computed using EOS BSK21 with a model with 
M = \AM Q and C = 0.1645). Evidently, the use of the 
late-inspiral cut-off frequency /contact calls for a formal- 
ism able to describe the phasing up to the merger (here, 
the EOB formalism and its accurate high PN expanded 
representation discussed in the previous section). 

In Eq. (53) we have included also a parameter (3 as- 
sociated to the spin-orbit interaction and a parameter a 
associated to the spin-spin one [36]. These parameters 
are equal to 

P = ^ (113X1 + 75^) L • a A + (A <-> B), (64) 

a = il (-247a A • a B + 721 L • k A L • a B ) , (65) 

where a a = $>a/{GM\) is the dimensionless spin param- 
eter of body A. Previous work [9, 31, 32] discussing data- 
analysis including the spin parameters j3 and er had incor- 
porated Baycsian priors a la [31] constraining the mag- 
nitudes of \/3\ and \<r\ to be smaller than 8.5 and 5.0 re- 
spectively, which are plausible theoretical upper limits on 
them. However, such values are very conservative bounds 
on p and a in view of observed binary pulsar systems (as 
already pointed out in Refs. [31, 36]). Indeed, recent 
estimates of the event-rate for BNS GW observations 
are mainly obtained from extrapolation of the currently 
observed binary pulsar systems. All the known binary 
pulsar systems have rather small observed spin parame- 
ters. Considering the fastest spinning pulsar observed in 
a BNS system, namely PSR J0737-3039A, whose spin pe- 
riod is 23 ms [37], we concluded from the calculations of 
moments of inertia by Bejger et al. [38] (who work with 
the EOSs: BPAL12, APR, SLy, BGN2H1 and GNH3) 
and by Morrison et al.[39] (who use FPS), that the ini- 
tial dimensionless spin parameter a is between approxi- 
mative^ 0.017 (for BPAL12) and 0.03 (for GNH3). This 
leads to an initial range for the corresponding parame- 
ter of order |/3| e [0.11; 0.196], while the 2PN-level 
spin-spin parameter a is at most of the order |<j| < 10~ 4 . 
Taking into account the slowing down of the spin until 
the moment of merger, we estimated that P at the time 
of the merger would be within the range [0.09; 0.17] so 
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that we decided to use the conservative upper limit of 
0.2 for /3. Hence, we studied the measurability of /? to- 
gether with the five other parameters {Ai,...,4; At} sub- 
mitted to a Gaussian Bayesian prior oc exp[— 1/2(/3/0.1) 2 ] 
constraining \j3\ to be smaller than 0.2 at the 95% confi- 
dence level. The result of the error estimates coming from 
such a constrained, six-parameter Fisher matrix formal- 
ism will be presented in Table II below, where they are 
compared to the result of a five-parameter Fisher ma- 
trix formalism where (3 is set to zero from the beginning. 
One sees from the numbers in Table II that such a con- 
strained six-parameter analysis leads to only a very slight 
increase of the error estimates. In view of this, in the fol- 
lowing we shall neglect (i.e. set to zero from the start) 
j3. Similarly, and a fortiori, in view of the very small 
upper bound quoted above on er, wc can also neglect the 
2PN level spin-spin parameter a. Let us emphasize that 
if, by contrast, one keeps the parameter /3 while using 
the very conservative prior j3 < 8.5, this leads to a very 
large increase of the error bars on M. and v 1 and a no- 
ticeable increase of the error bar on At- As it will be 
exemplified in Table II, the use of the very conservative 
prior constraining f3 < 8.5 instead of the "realistic" one 
leading to j3 < 0.2, increases the statistical measurement 
error on G/i2 by a factor which varies between 1.28 (for 
EOS BSK19) and 1.10 (for EOS GNH3). In addition, if 
one does not neglect the spin-spin parameter a (as we 
shall do here), or alternatively, does not put a realistic 
prior on it, but instead fits for it using a seven parame- 
ter Fisher matrix, constrained by the very conservative 
bound \a\ < 5.0, this leads to a further, substantial in- 
crease of the measurement errors. 



IV. MEASURABILITY OF TIDAL 
PARAMETERS: NUMERICAL RESULTS 

A. A sample of equations of state 

In this paper we consider a sample of EOSs taken from 
the literature. The sample is chosen to include EOS 
with a large range of variation in radius R, Love num- 
ber &2 and tidal parameter Gfi2- We consider eleven 
state-of-the-art EOSs. Seven among them, namely, MSI, 
MS2 [40], MPA1, AP3 [41], APR, SLy and FPS, have nor- 
mal matter content (npefi). One, namely GNH3 [42], also 
incorporates some mixture of hyperons, pion condensates 
and quarks. Finally, the three labels BSK19, BSK20 
and BSK21 refer to Skyrme-force-related energy density 
functionals (fitted to nuclear mass data) from which one 
can compute the EOS of cold neutron star matter [43]. 
Among these equations of state, seven of them (MSI, 
MS2, MPA1, AP3, SLy, FPS and GNH3) have been used 
in Ref. [9]. Table I lists, by order of decreasing radius 
(or increasing compactness) the main characteristics of 
Tohnan-Oppcnhcimcr-Volkoff neutron star models built 
from these EOS having mass 1.4M Q . These NS prop- 
erties were computed starting from the tabulated EOS, 



using Hermite polynomials interpolation [44], for all EOS 
but (MSI, MS2, MPA1, AP3), that we read instead from 
Table I of Ref. [9] . In our Table I, k^ refers to a fiducial 
1.4M Q + 1.4M Q binary system and / c = / contact refers to 
the contact frequency defined above (see Eq. (36)). Note 
that k,2 and Gfi2 decrease correlatively with the radius 
due to the dominant influence of the fifth power of the 
radius in kJ an d G^2 , and in spite of the nonmonotonic 
behavior of fe. The fifth root of Gfi2 defines a length 
scale which we can call the tidal radius of the NS. It is 
related to the radius R according to 
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(66) 



The values of the tidal radius for the 1.4Mq models 
listed in Table I vary between 8.8195 km (for MSI) and 
5.7297 km (for BSK19). The median value is around 
7 km. In the following we shall focus on a subsample 
of the EOS listed above, namely wc shall consider only 
GNH3, BSK21, BSK20, SLy, APR, FPS and BSK19, 
which span a plausible subrange of values of G/j.2 (note 
that we conservatively eliminate for instance the very stiff 
EOS MSI which yields an extremely large value of G^)- 



B. Measurability of Gfi2- equal-mass case 

In this section we focus on the measurability of tidal 
parameters in equal-mass BNS systems. We shall see 
below that, within the reasonable range of mass ratios 
expected from observational data, this equal-mass study 
is a sufficiently accurate indicator of the general case. 

Among the EOSs listed in the previous section, two 
of them (FPS and BSK19), which would have given the 
two smallest values of G^2, lead to maximum neutron 
star masses which are smaller than the recently reported 
value [25] M NS = (1.97 ± 0.04)M o . Because of this 
we shall first discuss the measurability of tidal param- 
eters within the restricted, observationally compatible, 
EOS subsample GNH3, BSK21, BSK20, SLy, APR. For 
each of these EOS we computed the 5x5 reduced Fisher 
matrix Fij, Eq. (47), corresponding to the parameters 
[Ai, . . . , A4; At] where the first four parameters refer to 
the binary system (see Eqs. (55)-(58)), while the tidal 
parameter At, defined in Eq. (59), reduces simply to 
G/i2 = G/j.2 = G/j.2 in the equal-mass case. The compu- 
tation of the Fisher matrix elements is performed by con- 
sidering that the GW signal is cut off above the (compact- 
ness dependent) contact frequency, Eq. (36), i.e. each 
integral is taken over the frequency window [/ m i n , /max], 
with / min = 10 Hz and / max = /contact- 

The diagonal elements of the inverse of the matrix F^ 
yield, according to Eq. (51), the SNR- normalized proba- 
ble (statistical) errors a\ i = p<r\ l on each parameter A^. 
Before discussing the measurability of the nontidal pa- 
rameters, let us start by considering the measurability of 
the tidal parameter At = G/i2- 
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TABLE I. Properties of a Ma = 1.4M© neutron star for the 11 EOSs considered in this paper. From left to right, the columns 
report: the name of the EOS; the radius of the star R; its compactness C; the value of the £ — 2 Love number fe; the value 
of the tidal parameters kJ and Gp,2\ the value of the bare contact frequency /bare m H z an d the corresponding dimensionless 
contact frequency /bare = /bare//o with /o = 56.56 Hz. 



EOS 


R [km] 


C 


k 2 


re 2 


G/x 2 [km 5 ] 


F [Hz] 


r 


MSI 


14.92 


0.1390 


0.1100 


264.9895 


53360.36 


1196.09 


21.15 


GNH3 


14.19 


0.1457 


0.0852 


162.099 


32641.6 


1284.06 


22.702 


MS2 


13.71 


0.1510 


0.0883 


140.6002 


28312.35 


1354.28 


23.94 


BSK21 


12.57 


0.1645 


0.0930 


96.4647 


19424.9 


1540.29 


27.23 


MPA1 


12.47 


0.1660 


0.0924 


91.6308 


18451.49 


1561.00 


27.60 


AP3 


12.09 


0.1710 


0.0858 


73.3528 


14770.89 


1632.06 


28.85 


BSK20 


11.75 


0.1760 


0.0810 


59.8628 


12054.4 


1704.83 


30.14 


SLy 


11.74 


0.1766 


0.0767 


55.8421 


11244.8 


1712.55 


30.2785 


APR 


11.37 


0.1819 


0.0768 


48.2159 


9709.13 


1790.2 


31.6514 


FPS 


10.85 


0.1907 


0.0662 


32.7966 


6604.17 


1922.7 


33.99 


BSK19 


10.75 


0.1924 


0.0647 


30.6662 


6175.14 


1948.14 


34.444 



A recent summary [1] of the expected event rate of BNS 
coalescences suggests that at the standard SNR detection 
threshold p = 8 a "realistic" estimate of the number of 
events per year detectable by the advanced LIGO- Virgo 
network is ~ 40. This means that at the SNR p = 16 one 
can reasonably expect to detect ~ 40/ (2 3 ) = 5 events per 
year. Considering such a SNR p = 16 we plot in Fig. 4, 
for each of the five EOS selected above, the following two 
curves: (i) the value of the tidal parameter Gp-i (in [km ]) 
as a function of the mass of each NS (thick, solid lines) ; 
and (ii) the corresponding value of the absolute statistical 
error oq^ (in [km ]). To guide the eye, a vertical line 
indicates the "canonical" mass value M = IAMq. If 
we first focus on this mass value, this figure shows that 
a single advanced LIGO or Virgo detector can measure 
Gp 2 for all considered EOS (with M max > 1.97M Q ) at 
a signal to noise ratio G/^/ctg/^ that varies between 1.4 
for the APR EOS up to 3.1 for the GNH3 one 5 . 

For mass values smaller than 1.4M Q the measurabil- 
ity of G/i2 is even better (larger ratio between Gp^ and 
ffcuJ, while for mass values larger than IAMq the mea- 
surability degrades. The intersection points in Fig. 4 be- 
tween solid and dashed curves corresponding to the same 
EOS mark the value the mass where Gpi is only measur- 
able at the "ler" (68% confidence) level, i.e., (Jg^ 2 = Gp2- 
For instance, for the APR EOS, still assuming a SNR 
p = 16, equal-mass BNS systems with individual NS 
masses larger than 1.52M© cannot allow one to extract 



5 These measurability ratios refer to observations by a single de- 
tector. Observing the same individual BNS event with a network 
of 3 LIGO- Virgo detectors will improve the measurability by a 
factor of order \/3 = 1.73, thereby leading to signal to noise ra- 
tios G/^2/o"g M2 varying between 2.4 for the APR EOS up to 5.4 
for the GNH3 one. 



GpL2 at a significant level. By contrast, in the case of 
BSK21 and GNH3 EOS one can extract tidal parame- 
ters for BNS systems up to individual masses larger than 
about 1.74M@. 

In summary, Fig. 4 shows that gravitational wave ob- 
servations from a single advanced detector are able to 
extract tidal parameters at a significant level, even for 
the soft EOSs that lead to the smallest values of Gp2- 
This conclusion strikingly contrasts with that of Hin- 
derer et al. [9] . We will discuss below the reasons behind 
this difference in conclusion. 

To complement the graphical representation of our re- 
sults in Fig. 4, we present in Table II numerical data 
referring not only to the 5x5, spinlcss, Fisher ma- 
trix calculation behind this figure, but to other calcu- 
lations. More precisely, this table gives SNR-normalized 
errors a\ t for all parameters of direct physical signifi- 
cance, namely A3 = M, A4 = v and At = Gp.2- [ Note 
that the numerical value of each d\ i formally gives the 
error corresponding to a unit SNR, p = 1. For larger 
values of p the error has to be divided by p.]. This table 
now considers the larger sample of EOS made by GNH3, 
BSK21, BSK20, SLy, APR, FPS and BSK19. For each 
one of these EOS we computed a 6 x 6 (or 5x5, see 
below) reduced Fisher matrix Fij, Eq. (47), correspond- 
ing now to the parameters [Ai, . . . , A4; /3; At]- Here, in 
addition to the first four binary-system parameters con- 
sidered above and of the tidal parameter At we also con- 
sider the spin-orbit parameter j3 (which will be treated 
with various different constraints, see below). We use the 
same frequency window ( [/ min , / max ], with / min = 10 Hz 
and / m ax = /contact) as above. We now consider the 
diagonal elements of the inverse of the matrix F\j for 
all the parameters of direct physical significance, namely 
A3 = M., A4 = v and At = Gp2, we list in Table II the 
corresponding SNR-normalized errors a\ i . 
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TABLE II. Measurability of the tidal parameter Gfi2 for a M = IAMq + IAMq neutron star binary obtained using the TaylorF2 
frequency-domain approximant to the phase, truncated at 2.5PN fractional accuracy for the tidal part and at 2PN accuracy in 
the point-mass part. From left to right the columns report: the name of the EOS; the value of the spin-orbit parameter and of 
the prior on it; the radius of the star; its compactness; the value of the £ — 2 tidal parameter Gfi2', the SNR- normalized relative 
errors on the chirp mass ctm/M and on the v parameter a v /u when cutting at bare contact frequency; the SNR- normalized 
absolute error, (Jg^.^, and relative error <7G/i 2 /(G7^2) on Gui- The last two columns refer to the absolute and relative errors on 
G\ii that are obtained by taking as cut-off frequency the conservative value 450Hz. 



EOS 


/? 




R [km] 


C 


G(J, 2 [km 5 ] 


0\nM 


<5"lny 


<7g P2 [km 5 ] 


<5"ln G/J 2 


*g?f [km 5 ] 


- 450Hz 
CT ln GM2 


GNH3 


1/3 


< +CXJ 


14.19 


0.1457 


32641.6 


0.00415853 


3.18959 


186 292 


5.70720 


1 476 380 


45.23 




1/3 


< 8.5 


14.19 


0.1457 


32641.6 


0.00405962 


3.09906 


182 612 


5.59447 


1 236 580 


37.8835 




1/3 


< 0.2 


14.19 


0.1457 


32641.6 


0.000447397 


0.122751 


165 714 


5.07679 


874 001 


26.7757 




P-- 


= 


14.19 


0.1457 


32641.6 


0.000450135 


0.117804 


165 652 


5.07487 


873 019 


26.7456 


BSK21 


1/3 


< +oo 


12.57 


0.1645 


19424.9 


0.003946 


2.98317 


158 080 


8.13801 


1 539 610 


79.2596 




1/3 


< 8.5 


12.57 


0.1645 


19424.9 


0.0038749 


2.91796 


155 190 


7.98922 


1 284 240 


66.1132 




1/3 


< 0.2 


12.57 


0.1645 


19424.9 


0.000434397 


0.115657 


133 108 


6.85246 


876 337 


45.1141 




P~- 


= 


12.57 


0.1645 


19424.9 


0.000436901 


0.110806 


133 046 


6.84928 


875 290 


45.0603 


BSK20 


1/3 


< +OC 


11.75 


0.1760 


12054.4 


0.00384331 


2.88426 


148 380 


12.3092 


1 575 360 


130.687 




\P 


< 8.5 


11.75 


0.1760 


12054.4 


0.00378349 


2.82927 


145 750 


12.0910 


1 311 380 


108.788 




1/3 


< 0.2 


11.75 


0.1760 


12054.4 


0.000428026 


0.112247 


118 815 


9.85656 


877 640 


72.8064 




P~- 


= 


11.75 


0.1760 


12054.4 


0.000430414 


0.107437 


118 751 


9.85125 


876 558 


72.7166 


SLy 


1/3 


< +oo 


11.74 


0.1766 


11244.8 


0.00383898 


2.8801 


148 911 


13.2426 


1 579 310 


140.448 




1/3 


< 8.5 


11.74 


0.1760 


11244.8 


0.00377961 


2.82552 


146 254 


13.0064 


1 314 390 


116.888 




1/3 


< 0.2 


11.74 


0.1760 


11244.8 


0.000427755 


0.112104 


118 271 


10.5179 


877 784 


78.0612 




/3 = 


= 


11.74 


0.1760 


11244.8 


0.000430139 


0.107295 


118 206 


10.5121 


876 697 


77.9646 


APR 


1/3 


< +OC 


11.37 


0.1819 


9709.13 


0.00379747 


2.84028 


142 857 


14.7136 


1 586 810 


163.434 




1/3 


< 8.5 


11.37 


0.1819 


9709.13 


0.00374226 


2.78947 


140 408 


14.4615 


1 320 100 


135.964 




1/3 


< 0.2 


11.37 


0.1819 


9709.13 


0.000425161 


0.110728 


112 643 


11.6018 


878 055 


90.436 




/3 = 


= 


11.37 


0.1819 


9709.13 


0.000427498 


0.105935 


112 580 


11.5953 


876 961 


90.3233 


FPS 


1/3 


< +oo 


10.85 


0.1907 


6604.17 


0.00373437 


2.77992 


135 473 


20.5133 


1 602 010 


242.575 




1/3 


< 8.5 


10.85 


0.1907 


6604.17 


0.00368509 


2.73448 


133 267 


20.1792 


1 331 690 


201.644 




1/3 


< 0.2 


10.85 


0.1907 


6604.17 


0.000421197 


0.108641 


104 424 


15.8118 


878 605 


133.038 




/3 = 


= 


10.85 


0.1907 


6604.17 


0.000423462 


0.103871 


104 362 


15.8025 


877 496 


132.87 


BSK19 


1/3 


< +OC 


10.75 


0.1924 


6175.14 


0.00372323 


2.76928 


134 005 


21.7007 


1 604 110 


259.769 




1/3 


< 8.5 


10.75 


0.1924 


6175.14 


0.00367495 


2.72475 


131 846 


21.3511 


1 333 300 


215.914 




1/3 


< 0.2 


10.75 


0.1924 


6175.14 


0.000420494 


0.108273 


102 998 


16.6795 


878 681 


142.293 




/3 = 


= 


10.75 


0.1924 


6175.14 


0.000422746 


0.103507 


102 937 


16.6696 


877 570 


142.113 



For each EOS, the results are displayed along four 
rows. On each row, the first four columns give: (i) infor- 
mation about the treatment of the spin-orbit parameter 
/3; (ii) the value of the neutron star radius (in km); (iii) 
the value of the compactness; (iv) the value of the tidal 
parameter G/j,2- The following four columns give: (v) 
the fractional, SNR normalized, error on the chirp mass 
&\ n M = &mIM.; (vi) the fractional SNR normalized, er- 
ror on the symmetric mass ratio, a\ uv = o v jv\ (vii) the 
absolute, SNR normalized error &g^ 2 on G[ii (in [km ]); 
and finally (viii) the fractional, SNR normalized error 
ai n Gfj,2 — <5'GM2/(G'/ i 2) on G/j.2- Concerning the treat- 
ment of the spin-orbit parameter, the first row, labelled 
with |/3 1 < +oo refers to a 6 x 6 Fisher matrix analysis 



where /3 is included as a sixth unconstrained parameter. 
The second row, |/3| < 8.5, refers to a a 6 x 6 Fisher ma- 
trix analysis where /3 is constrained by adding a Gaussian 
prior proportional to exp[— 1/2(2/3/8. 5) 2 ]. Similarly, the 
third row corresponds to a more constraining prior pro- 
portional to cxp[-l/2(2/3/0.2) 2 ]. Finally, the fourth row 
corresponds to a 5 x 5 Fisher matrix analysis where /3 is 
set to zero from the beginning without being fitted for, 
which was used to obtain the data displayed in Fig. 4. As 
already mentioned above, the results for the strong prior 
|/3| < 0.2 (3rd row) are nearly indistinguishable from the 
results of the 5x5 Fisher matrix analysis (4th row) . This 
justifies our use of the 5x5 Fisher matrix results in Fig. 4 
above. By contrast, we see that the results corresponding 



15 



a. 
Q 



. . , ' ' ' ; 


\_ ^ 


; ^\ p=w ; 


BSK21 X. 


"— — — . ^^^ N. bsfe21 y. 


'. "~~-~~~~~^^^ N. gnh3 _ _^~- "'bsk20 ^ 


^ 3^ """*-"-' -*""**" 


- BSK20 ^^^^^-lI^C " --*■"*" " 






SLY^*~~--^C^---r^^^_, - " " " , - ' ' = = = S ' ' ^^---^ 


~~ ^ -— - s si- s ;l p r ^v^"**-^-« 11 ^ iii ^ 






,,-:;^ = = = ~ ^ ====: ^^^r^~~^. N. 




^=^ 



1.6 1.7 1.8 1.9 

Mass (M e ) 



FIG. 4. Measurability of the tidal polarizability parameter 
G/i2 (in units of km 5 ) as a function of the neutron star mass 
for a sample of realistic EOS from Table I. This plot refers to 
the observation (at the SNR level p = 16) of the gravitational 
wave signal from an equal-mass BNS merger as seen by a sin- 
gle advanced LIGO detector. The solid lines represent the 
values of Gp2 as a function of the NS mass, while the dashed 
lines represent the 1 a (68% confidence level) expected sta- 
tistical errors. The vertical line marks the canonical NS mass 
IAMq. Note that over a wide range of masses each solid 
line lies comfortably above the corresponding measurability 
threshold, therefore indicating that the advanced LIGO- Virgo 
detector network can significantly measure Gp2- 



either to the conservative prior \j3\ < 8.5 (second row) or 
the lack of any prior (first row) are close to each other 
but differ from the strongly /3-constrained results by very 
significant factors. To be precise, the measurability of the 
chirp mass is worsened by a factor larger than seven; that 
of the symmetric mass ratio is worsened by a factor of 
order 30!; finally, that of G/i2 is only worsened by about 
20%. These results arc linked to the different origins of 
the effective signals contributing to the measurability of 
the various parameters displayed in Fig. 3. 

We can roughly summarize the results for the measur- 
ability of the nontidal parameters (in the strongly con- 
strained (3 cases) in the following way: 



excellent (and has not been very significantly worsened 
by the inclusion of the tidal term, as shown by compar- 
ing to the results of Refs. [31, 32]). By contrast, the 
fractional precision on v has been significantly worsened 
(by a factor of order 1.7) compared to Refs. [31, 32] when 
fitting for an extra tidal parameter 6 . This worsening in 
the measurability of v might make it difficult to distin- 
guish stars with a mass ratio between 0.75 and 1. For 
instance, if wc considered a BNS with Ma = 1.2M Q , 
M B = 1.6M (i.e., M A /M B = 0.75) its symmetric mass 
ratio is v w 0.2449, so that 1 - Av = 0.0204, correspond- 
ing to a fractional 5v/v w 0.02. Comparing this with the 
measurement error in v for p = 8, Eq. (68), this is only a 
2cr-level deviation. Actually, this problem may be cured 
by doing two separate analyses of the GW data, one using 
inspiral data only up to a cut-off frequency small enough 
to be able to neglect tidal effects (without trying to fit 
for tidal parameters), which will probably give a better 
estimate of the mass ratio. And a separate analysis of 
the data up to (and possibly beyond) the merger aimed 
at extracting EOS-dcpendent information. 

The last two columns of the table exhibit the SNR- 
normalizcd absolute and relative errors on Gp2 in the 
case where one uses as upper frequency cut-off / max = 
450 Hz as done in Ref. [5, 9]. The use of such a lower 
cut-off leads to a dramatic worsening (by a factor ~ 7) 
of the measurability of G\i2 (the origin of this worsening 
is illustrated in Fig. 2, which includes a line at 450 Hz). 

On the other hand, Hindercr ct al. [9] computed a 
SNR- normalized uncertainty on Gpt2 for the 1.4M Q + 
1AM Q system equal to <7§)? 2 derer = 35 x 19.3 x 
0.6674310 4 km 5 = 450.84 x 10 4 km 5 (see second row of 
their Table II which corresponds 7 to a SNR p = 35). 
Considering for example the SLy EOS, this is a factor 
38 larger than the corresponding result in Table II for 
our preferred 5-parameter analysis. This large factor 
can be viewed as originating from the product of sev- 
eral subfactors: (i) a factor of order (/ c /450 Hz) 2 - 2 = 
(1704/450) 22 w 18.7 due (according to Eq. (23) of 
Ref. [9]) to their use of a cut-off at 450 Hz; (ii) a fac- 
tor ~ 1.24 due their use of a conservative prior (8.5) on 
/3; iii) a supplementary factor coming from the fact they 
also fit for the 2PN spin-spin parameter a (with a con- 
servative prior), thereby working with seven correlated 
parameters. 



M 



and 



v 



4.3 x 10" 
P 



0.11 
9 



(67) 



(68) 



For instance, when p = 10 this means that the chirp mass 
is measured to a fractional precision of 4 x 10 -5 , while 
the symmetric mass ratio is measured at a fractional pre- 
cision of 0.01. As usual, the fractional precision on A4 is 



6 Note that when one is fitting for the spin parameter /3, the frac- 
tional precision of v becomes dramatically worsened, down to 
the level a\ nv ~ 2.8. In the case of EOSs GNH3 and BSK21 this 
renders the fractional accuracy on v comparable to the fractional 
accuracy on G/j.2 ■ In such a case there can be a large difference 
in the measurability of At, Eq. (59) versus X' T , Eq. (61), espe- 
cially in view of the correspondingly large correlation between 
G/i2 and v. 

Wc could not reconcile the statement in Ref. [9] that they con- 
sider a source at a distance of 100 Mpc, with an amplitude av- 
eraged over sky position and relative inclination, with the SNR 
35 quoted in their Table II, which, according to Abadie et al. [1] 
seems to correspond to an optimally oriented source at 100 Mpc. 
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TABLE III. Measurability of the tidal parameter pa, Eq. (59), for a M = 1.1529M Q +1.6470M Q neutron star binary (M A /M B = 
0.7) obtained using the TaylorF2 frequency-domain approximant at 2.5PN fractional accuracy in the tidal part of the phase 
and at 2PN accuracy in the point-mass part of the phase. From left to right the columns report: the name of the EOS; the 
compactnesses of the stars; the values of the £ = 2 tidal polarizability parameters G\ia and G\xb of the stars; the value of the 
tidal parameter Gp,2', the bare contact frequency /b arc in Hz; the SNR-normalized relative errors on the chirp mass &m/M 
and on the v parameter g v Jv\ the absolute (in km 5 ) and relative errors on Gp.2- For each EOS, the second row recalls the 
corresponding results for equal masses. 



EOS 


Ca 


C B 


G/i2 4 [km 5 ] 


G^[km 5 ] 


GfL 2 [km 3 ] 


/bare [Hz] 


<5"ln M 


O'bm 


<7Gp2 [km 5 ] 


<5"ln Gp, 2 


GNH3 


0.118 


0.178 


41115.1 


20711.1 


36178.1 


1301.97 


0.000429438 


0.115405 


158 647 


4.38517 


GNH3 


0.1457 


0.1457 


32641.6 


32641.6 


32641.6 


1284.06 


0.000450135 


0.117804 


165 652 


5.07487 


BSK21 


0.1361 


0.1938 


22049.1 


15623.7 


21034.5 


1547.11 


0.000417681 


0.108921 


129 240 


6.1442 


BSK21 


0.1645 


0.1645 


19424.9 


19424.9 


19424.9 


1540.29 


0.000436901 


0.110806 


133 046 


6.84928 


BSK20 


0.1446 


0.2091 


14638.3 


8891.9 


13429.4 


1713.76 


0.00041153 


0.105588 


115 330 


8.5879 


BSK20 


0.1760 


0.1760 


12054.4 


12054.4 


12054.4 


1704.83 


0.000430414 


0.107437 


118 751 


9.85125 


SLy 


0.144 


0.2116 


14347.3 


7696.2 


12794.1 


1723.62 


0.000411201 


0.105412 


114 634. 


8.95991 


SLy 


0.1760 


0.1760 


11244.8 


11244.8 


11244.8 


1712.55 


0.000430139 


0.107295 


118 206 


10.5121 


APR 


0.14934 


0.2157 


11874.0 


7144.9 


10868.9 


1797.14 


0.000408861 


0.104155 


109 549 


10.0792 


APR 


0.1819 


0.1819 


9709.13 


9709.13 


9709.13 


1790.2 


0.000427498 


0.105935 


112 580 


11.5953 


FPS 


0.154 


0.2345 


9443.33 


3433.1 


7830.8 


1956.63 


0.000404363 


0.101758 


100 125 


12.7862 


FPS 


0.1907 


0.1907 


6604.17 


6604.17 


6604.17 


1922.7 


0.000423462 


0.103871 


104 362 


15.8025 


BSK19 


0.1553 


0.2352 


8866.4 


3386.96 


7411.5 


1973.3 


0.000403933 


0.101529 


99 261 


13.3929 


BSK19 


0.1924 


0.1924 


6175.14 


6175.14 


6175.14 


1948.14 


0.000422746 


0.103507 


102 937 


16.6696 



C. Measurability of G/12: unequal-mass case 

Let us now consider the measurability of tidal parame- 
ters in unequal-mass BNS systems. Following Refs. [9, 45] 
we focus on comparing the measurability of G/i2 in a sys- 
tem with a large, but plausible, mass ratio Ma/Mb = 0.7 
(corresponding to v = 0.2422) to an equal-mass system. 
Taking the total mass of the system to be the canon- 
ical M = 2.8M Q , the mass ratio we chose determines 
M A = 1.1529M Q and M B = 1.6470M©. Here we use the 
5-paramctcr Fisher- matrix analysis with /3 = 0. For the 
same sample of EOS as in Table II, Table III lists the 
individual compactnesses, the tidal parameters Gfi 2 ' , 
their combination G/12 and the SNR-normalized abso- 
lute uncertainty on G/22 , ^Gp, 2 as weu as the relative one 
^inGp.2- For improved readability of the table, for each 
EOS we also include the equal-mass result of Table II in 
a second row. 

As the numbers in this table show, even the large mass 
ratio 0.7 that we consider does not influence much the 
measurability of tidal parameters. In all cases it im- 
proves both the absolute and the fractional measurabil- 
ity of G/2 2 by about 10 - 15% (for the EOSs that lead to 
Mmax > 1.97M , i.e. excluding FPS and BSK19). 

Note that the computations behind the results in Ta- 
ble II-III have assumed a cut-off frequency which was a 
function of the individual compactnesses of the two stars. 
In an actual GW data analysis situation we won't have an 
a priori knowledge of these compactnesses and therefore 



we will need a way to internally fix the value of the fre- 
quency up to which the EOB template can be considered 
as a reliable description of the observed GW signal. One 
can think of several ways in which this could be done. 

A first way is to use the fact that, for a given EOS, 
the contact frequency J contact ; or the contact frequency 
parameter x C ontact given by Eq. (36) is symmetric under 
the A <->• B exchange, and therefore it can (in principle) 
be considered as a function of Gp,2, M, and v (which are 
all A f-> -B-symmetric functions). Moreover, in view of 
the strong, approximately universal, dependence of G[i 2 
on C A (Gfi A ex k£(C A )CA 5 , with h£(C A ) « A - BC A , 
as per Eq. (116) of [6], using for (A,B) values appro- 
priate for an "average" EOS) the function determining 
^contact j n t erms f G/22, M and v can be considered as 
approximately universal and known. A second way is, 
separately from a tidal-parameter-fitting data analysis 
of the inspiral signal, to use the full GW data includ- 
ing the post-merger signal to extract information both 
about the frequency of merger and the post-merger dy- 
namics, so as to have some independent handle on the 
EOS. Indeed, recent numerical results [46-49] on BNS 
merger have shown that the GW signal contains definite 
imprints both of the merger and post-merger dynamics. 
For instance, on Fig. 2 of Ref. [49] both the frequency 
marking the end of inspiral (corresponding to J contact m 
our EOB setup) , and the characteristic frequency of post- 
merger oscillations stand out above the advanced LIGO 
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TABLE IV. Correlations, as given by Eq. (69), between 
the dynamically relevant parameters, {M,v,Gfi2} of the 
1.4Mq + IAMq binaries of Table II as obtained from the 
5x5 Fisher matrix analysis (the spin parameters, (/3, a), are 
set to zero from the start). From left to right the columns re- 
port: the name of the EOS, the cut-off frequency (either con- 
tact frequency or 450 Hz), the correlation between (Gfj,2,M), 
the correlation between (Gfi2, v) and the correlation between 
(M,v). The correlations between Gfi2 and (M,v) decrease 
when the cutoff frequency is increased. 



EOS 


/ max [Hz] 


°35 


r ,(v,Gfj. 2 ) 
°45 


MM,*) 

°34 


GNH3 


1284.06 


0.568314 


0.755652 


0.911907 


GNH3 


450 


0.693511 


0.859695 


0.93168 


BSK21 


1540.29 


0.551751 


0.73937 


0.909245 


BSK21 


450 


0.694515 


0.860478 


0.93168 


BSK20 


1704.83 


0.543526 


0.730951 


0.907868 


BSK20 


450 


0.695072 


0.860911 


0.93168 


SLy 


1712.55 


0.543398 


0.730791 


0.907809 


SLy 


450 


0.695133 


0.860959 


0.93168 


APR 


1790.2 


0.539288 


0.726584 


0.907232 


APR 


450 


0.695249 


0.861049 


0.93168 


FPS 


1922.7 


0.539288 


0.726584 


0.907232 


FPS 


450 


0.695249 


0.861049 


0.93168 


BSK19 


1948.14 


0.532037 


0.719014 


0.906174 


BSK19 


450 


0.695515 


0.861256 


0.93168 



D. Correlations 

To complement the results about the measurability of 
G/Z2 (and G\±2) given in the previous two sections, let us 
discuss the issue of the correlations among the various 
parameters and their influence on the measurability of 
G/i2- Usually, correlations are measured via the nondi- 
agonal terms of the covariance matrix, that is by 



-iyj 



Cj. 



(F- 1 ) 



rp-iyifp-iyj 



(69) 



which arc numbers that vary between —1 and +1. We 
focus on canonical 1.4M Q + 1AM@ equal-mass binaries 
built from our sample of EOSs. The values of the corre- 
lations (for the dynamical parameters {M.,v, G/^2}) Cij 
for the 5x5 Fisher matrix analysis are listed in Table IV. 
For each EOS, we list the values of Cij when one takes as 
cutoff frequency the contact frequency (top row) as well 
as their values when taking 450 Hz as cutoff frequency 
(bottom row). Note first that G/i2 is only (especially 
when using the contact frequency as cutoff) moderately 
correlated to M. and v. by contrast to the (M.,v) correla- 
tion, the (M,G/j,2) and (v, G/X2) correlations are always 
comfortably smaller than 0.9. For a given cut-off fre- 
quency, the values of the correlations Cij decrease when 
the compactness of the model increases. This decrease 



TABLE V. The quantities G t , as defined by Eq. (70), for the 
physically relevant parameters, {M,v, G/^2} of the IAMq + 
1.4Mq binaries of Table II. These values are obtained from 
the 5x5 Fisher matrix analysis, where the spin parameters, 
(/3,a), are set to zero from the start. From left to right the 
columns report: the name of the EOS, the cut-off frequency 
(either contact frequency or 450 Hz), and the Gi's. The value 
of d decreases when the the cutoff frequency is increased. 



EOS 


/ max [Hz] 


q(M) 


Ur 4 


G (G M2 ) 


GNH3 


1284.06 


3.80973 


12.7126 


4.94471 


GNH3 


450 


5.04097 


21.973 


10.9479 


BSK21 


1540.29 


3.69651 


11.959 


4.44276 


BSK21 


450 


5.04099 


21.982 


10.9764 


BSK20 


1704.83 


3.6412 


11.5971 


4.20174 


BSK20 


450 


5.041 


21.987 


10.9923 


SLy 


1712.55 


3.63886 


11.5821 


4.19293 


SLy 


450 


5.04101 


21.9876 


10.994 


APR 


1790.2 


3.57594 


11.1741 


3.91152 


APR 


450 


5.04101 


21.9886 


10.9973 


FPS 


1922.7 


3.58202 


11.2133 


3.9388 


FPS 


450 


5.04101 


21.9908 


11.004 


BSK19 


1948.14 


3.57594 


11.1741 


3.91152 


BSK19 


450 


5.04102 


21.9911 


11.0049 



is mild. To be precise, considering the variability be- 
tween GNH3 and BSK19, we have that Gg 5 ' varies 
by 6.4%, C { £ G ^ ] by 4.8% and C^' v) by 0.6%. On the 
other hand, the values of the correlation increase when 
the cutoff frequency is decreased from the contact fre- 
quency to 450 Hz. This is expected since up to 450 Hz 
the tidal part of the phasing is quite weak and thus 
rather difficult to disentangle from the nontidal signal. 
Note that G/12 is only moderately correlated to M. and 
v. by contrast to the (M.,v) correlation, the (M,G/j,2) 
and (v, G/X2) correlations are always comfortably smaller 
than 0.9. 

It is also useful to look at the quantity (for each pa- 
rameter A,) 



G i = iJ(F-i)«F ii . 



(70) 



For each EOS wc list the values of Gi (for i = 3,4,5, 
i.e. Ai, v and G/12) in Table V. The quantity Gi mea- 
sures the global correlation coefficient, say cf ° a , of the 
parameters Xi with respect to all other parameters Xj, 
i ^ i, via 

1 



G,= 



(71) 



1 - 



c global\ ' 



Here. c : 



global 



is the larger possible correlation between 



Xi and a linear combination of the other parameters A, , 
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j y^ i. Let us discuss in more detail the meaning of the 
quantities G^. We recall that it is convenient to interpret 
the measurability of the various parameters Xi entering 
the phasing signal &(f; Aj) in terms of geometrical con- 
cepts related to the scalar product (49) (which is the SPA 
version of the Wiener scalar product). When considering 
small variations SXi of the parameters, to each parame- 
ter Xi is associated the vector &i = dx^ so that the in- 
finitesimal signal associated to a simultaneous variation 
of all the parameters is the following linear combination 
of individual vectors: S<P = ^\ SX^i. The geometrical 
transcription of the fact that the measurement of a par- 
ticular parameter Xi is correlated to the measurement of 
the other parameters Aj, j ^ i, is that the signal vec- 
tor $i associated to SXi is not orthogonal to the other 
signal vectors <Pj. [Remember that the Fisher matrix is 
the matrix of scalar products Fy = (<Z^|#j).] The global 
correlation between Xi and all the other Aj's, j =/= i, is 
then measured by the "inclination angle" 9i between the 
vector <Pi and the hyperplane Hi spanned by the remain- 
ing vectors \Pj, j ^ i. The angle 9i is defined so that 
it vanishes when Wi lies within the hyperplane Hi, and 
equals ir/2 when ^i is orthogonal to the hyperplane Hi. 
Let us now decompose the vector >?i in two orthogonal 
vectors: (i) its projection <pj orthogonal to Hi, and (ii) 

its projection tyf parallel to Hi. It is then easy to see 
that the definition of Gi given by Eq. (70) implies 



1 






(72) 



where \v\ denotes the (Euclidean) length of the vector 
v in signal space. Note also that the global correla- 
tion coefficient cf ° defined above is simply equal to 

global /i 

c\ =cos0i. 

Let us also note the following formulas yielding the 
SNR-normalized (absolute and fractional) error(s) on the 
parameter Xi 



<T\i = P°\, 



°ln A, 



X l 



1 

1^1 


G, 

m 1 


1 


G % 



M& 



Ai|#i 



(73) 



(74) 



In particular, if we apply the last formula to the tidal 
parameter Xt = G/j,2 we see that, given a certain SNR 
p, the two factors that determine the measurability of 



At are Gt = G 



(G M2 ) 



(which measures the adverse ef- 



fect of correlations with the other parameters, and which 
should be as small as possible) , and the Euclidean length 
of the full tidal signal |Ar^r|- In other words, the "use- 
ful" part of the explicit frequency-domain tidal signal 
At^t(/) pictured in Fig. 1 (which reaches about ten ra- 
dians at contact) is reduced by two factors: a first factor 
coming from the overlap between the SNR measure "/(f) 
and the tidal signal &r(f) ~ / 5 ^ 3 (which enters the in- 
tegral \\ T V T \ = J J df-y(f)(X T ^T(f)) 2 )), and a second 



factor 1/Gt, due to the correlations (which retain only 
the part of the tidal vector which is orthogonal to all the 
other signal vectors). This motivates us to define the 
useful number of radians (in a rms sense) contained in 
the tidal dephasing signal as 



^ sc 



fui 



AtI^tI 



X t \^t\ 
Gt 



1 



Cln G^2 



(75) 



In view of the results reported in Tables II and III a me- 
dian estimate for the useful tidal dephasing is i^ soful ~ 
0.1 rads. This is a factor ~ 100 smaller than the dephas- 
ing at contact. This reduction factor can be seen as the 
product of a factor Gt ~ 5 due to correlations, and a 
factor ~ 20 coming from the fact that the tidal signal is 
strongest during late inspiral, when the SNR curve f-f(f) 
is much below its maximum (see Fig. 3). Note also that 
if one uses 450 Hz as cut-off frequency the useful tidal 
dephasing is drastically reduced (roughly by a factor 7): 
e.g. for the BSK21 EOS which led to a rather comfort- 
ably measurable tidal signal 3^BSK2i ~ 1/7 when con- 
sidered up to contact, one has only ^bsk^l ~ 1/45- 
This loss in measurability by a factor ~ 6.6 is due both 
to a higher global correlation (Gt increasing from 4.44 
to 11.0) and to a smaller signal at 450 Hz versus f c . As 
in the case of the correlations Gy , the main message of 
Table V is that (especially when considering as cutoff 
frequency the contact frequency) the global correlation 
of Gfi2 with respect to all other parameters is moderate 
and comparable to that of M.. By contrast, v is more 
strongly correlated to the other parameters. 



E. Coherent data analysis of tidal parameters 

Until now we have been discussing the measurability of 
tidal parameters from the GW signal emitted by a single, 
particular BNS merger event (eventually simultaneously 
observed by 3 separate detectors). We wish now to in- 
troduce a new way of extracting EOS-dcpcndcnt infor- 
mation by a "coherent" data analysis of the GW signals 
emitted by many separate BNS merger events, say the 
expected ~ 40 BNS mergers observable in one year by 
one advanced LIGO (or Virgo) detector. [Evidently, the 
method can also be extended to a coherent analysis of 
the data coming from the full network of LIGO- Virgo 
detectors] . 

This method is based on the following preliminary re- 
mark. As exemplified on our Fig. 4 above (as well as 
in Fig. 2 of Ref. [9]), the tidal parameter G/i^ is, for a 
given EOS, a function of the mass Ma of the considered 
neutron star which can be well represented by a linear 
function in the range of expected neutron-star masses, 
1.2M Q - 1.9M Q , say 



Gi4(M A ) = a EO s + o-eosMa- 



(76) 



The crucial point here is that the coefficients (oeoSj ^eos) 
depend only on the EOS, but not anymore on the spe- 
cific neutron star mass. In the following, we shall use the 



19 



symbols (oeoSi ^eos) to denote the two unknown param- 
eters corresponding to the actual EOS chosen by Nature. 
Moreover, the system tidal parameter GJI2, Eq. (59), 
that enters the inspiral signal of an individual system 
(Ma, Mb), becomes, when using Eq. (76) 



1 /6 

G> 2 = t^oeos 11 

13 V v 



13 aEOS 



i-» 



2&eosM 

H 7T&EOS 



A 



3/5' 



(77) 



Let N denote the number of BNS merger events observed 
during a certain period T (e.g, 1 yr). We introduce an 
index / = 1 , . . . , N labelling each BNS system, and the 
corresponding merger event, within this collection of ob- 
served GW signals. We now discuss a data analysis pro- 
cedure for the combined event consisting of this collec- 
tion of N individual GW signals. This "grand signal" 
depends on a collection of parameters: (A 7 , oeoSi ^eos)- 
Here the A^'s vary from BNS system to BNS system (and 
include, besides the parameters Ai,...^ considered above, 
also an amplitude parameter), while the two EOS pa- 
rameters, (oeoSj &eos) are common to the whole col- 
lection of events. We now envisage a grand fit of the 
whole collection of parameters (A„, oeoS: ^eos) to the 
ensemble of N GW signals. The signals in this ensem- 

I 



ble have clearly statistically independent noise contri- 
butions. Let us then consider the Bayesian probabil- 
ity distribution function for the values of the parameters 
(A^, aEOB, &eos)j given a grand strain signal {s 7 (i)}. It 
can be written as exp (— x 2 /2) where [31] 



N 



X 2 (A^oeos^eos; s 1 , priors) = ^X/ pr 



1=1 



N 



2j(/i (A 7 , a E os,^EOs) ~ s I ,h(X I a ,a E os,bEOs) - s 1 ), 

(78) 



1=1 



with (...,...) denoting as above the single-observation 
Wiener scalar product, and X/ pi ior indicating the 
logarithms of eventual priors on some parameters 
(c-g-,X/ P rior = (2/3///3 m ax) 2 as above). For simplicity, 
let us use only strong priors, that are equivalent to elim- 
inating some parameters (e.g., /? max — > equivalent to 
setting Pi to zero). In the high SNR approximation, and 
after having marginalized over the N amplitude param- 
eters Ai (treated in the Gaussian approximation, and as 
approximately independent of the other parameters), we 
can approximately reexpress \ 2 in terms of phase differ- 
ences, using the renormalized scalar product (49): 



N 

E< 

1=1 



(Xi,a E os, &eos; s 1 ) - E ^ ( *"(/; A i, oeos, &eos) - *,/ (J)Mf\ x i, oeos, b 



EOS ) 



*Af) 



(79) 



The scalar product (49) a priori depends on the index 
I through the choice of the cut-off frequency / max = 
f c (M 1 ,C A ,X A ). As a first approximation for under- 
standing how using such a \ 2 improves the measure- 
ment of (aEOS,^EOs)j let us however consider that one 
uses some a priori fixed cut-off frequency. The theoreti- 
cal phase & (/; A£,a E os,&EOs) depends on (a E os,&EOs) 
only through a term of the form A^tidai(/; A 7 ), where 
A^ = G/2| is linear in a E os an d &eoSj see Eq- (77). If 
wc then project #tidai(/) into its projection &£ dI (f), 
which is orthogonal to the nontidal signals W^ 1 (f) — 
d\i l P (f;\ I a ), [with respect to the Euclidean metric 
(. . . I . . . )], we find that the part of \ 2 which depends 
on (dEOS, ^eos) is quadratic in them and of the form 



N 



X 2 (asos, b E os) -J2p 2 i (*tid j|#Hdl) 



1=1 



13VAT" 11,aEOS 



2(A|) 3 / 5 



&eos ~c(s ,A a ) 



(80) 



Here (as discussed in the previous section) the fac- 



I 

tor (^tidjl^tid/) takes into account the correlation of 
(oeoSj ^eos) (via A^) with the nontidal parameters. Fi- 
nally, the latter formula defines, for each confidence level, 
an error ellipse in the (aEOS, &eos) plane. The size of the 
minor axis of the ellipse (associated to some best deter- 
mined Ay-like combination of a E os an d &eos) w iH essen- 
tially be determined by the following effective squared 
SNR, given by the sum of all individual SNR's, i.e. 
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On the other hand the size of the major axis of the er- 
ror ellipse will crucially depend on the dispersion of the 
distribution of A3 = A4 and A4 = v around their median 
values. If such a dispersion is large enough, one will be 
able to measure both oeos and &eos with a precision 
still mainly determined by the effective SNR (81). If, 
on the contrary, all observed BNS systems happen to be 
close to, say, the canonical (1.4M Q , 1.4M Q ) system, the 
error ellipse will be very elongated in one direction, while 
the other direction (minor axis) will lead to the measure- 
ment of one Ay-likc combination of oeos and £»eos- ^ n a H 
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cases, the crucial quantity determining the measurability 
of some Ay-like quantity (or quantities) is the effective 
squared SNR, Eq. (81). One can then estimate the sum 
giving the effective squared SNR Eq. (81) by approxi- 
mating it by an integral over the ball of space around 
the Earth containing all the source events up to the min- 
imum threshold of dctectability of a BNS, that we shall 
take as p m - m = 8. This leads to the following estimate of 
the coherent SNR 



TABLE VI. Sensitivity of cti 



nG^2 



when varying a 



(2) 



/^coherent 



(82) 



where the fact that the p m in is augmented even beyond 
the naively expected \N factor is due to the existence of 
closer events having a larger SNR. If we apply this result 
to an expected "realistic" number of events, N ~ 40 dur- 
ing one year of observation, we conclude that the effective 
SNR for such a coherent analysis is ~ \/3 X 40 x8~ 88. 
On the other hand, there is clearly a price to pay for 
such an remarkable increase in SNR. Indeed, as sketched 
above, the two parameters (oeos, &eos) will now be cor- 
related among themselves, which will degrade their indi- 
vidual measurability. We leave to future work a detailed 
exploration of the performance of such a coherent anal- 
ysis (taking into account all the correlations that have 
been neglected in the sketchy treatment above), espe- 
cially for the measurability of (oeos, ^eos)- We however 
expect that it will significantly improve the single detec- 
tor measurability computed above, thereby strongly re- 
inforcing our conclusion that the advanced LIGO- Virgo 
network can extract EOS information from the late in- 
spiral BNS signal. 
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F. Sensitivity to a 



In the analysis of this paper we have used the recently 
analytically computed value of the 2PN tidal amplifica- 
tion parameter, Eq. (A10), leading to a 2 ' = 85/14 w 6 
in the equal-mass case. However, the comparison be- 
tween the waveform prediction coming from the EOB 
tidal model and recent, state-of-the-art BNS numerical 
simulations [22, 23] has suggested that the effective value 

(2) 

of a 2 might be, in the equal-mass case, of order 40 or 
even larger. In addition, analytical arguments have been 
advanced in Ref. [24] suggesting that a 2 might be fur- 
ther amplified by higher PN effects, possibly by a factor 
of order 2. In view of this uncertainty on the influence 
of higher relativistic corrections to the tidal interaction 
energy, we have explored the effect on the measurability 

(2) 

of Gp 2 of changing the value of a 2 . Our results arc dis- 
played in Tabic VI, which for a large sample of EOS lists 
the value of the SNR- normalized fractional error a\ n g^ 2 
on Gp 2 for canonical 1.4M© + 1.4M systems. As we 
see from the numbers in the table, if it turns out that 
the effective value of a 2 ' is closer to 40 than to 6, this 
will improve the measurability of Gfi 2 by more than 20%. 
Therefore, all our conclusions above should be considered 
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MSI 


3.5010 


3.3733 


2.7975 


GNH3 


5.2783 


5.0749 


4.1688 


MS2 


5.7074 


5.4774 


4.4629 


BSK21 


7.1710 


6.8493 


5.4649 


MPA1 


7.4385 


7.1011 


5.6529 


AP3 


8.8549 


8.4382 


6.6663 


BSK20 


10.3566 


9.8513 


7.7224 


SLy 


11.0534 


10.5121 


8.2342 


APR 


12.2162 


11.5953 


9.0079 


FPS 


16.7035 


15.8025 


12.1087 


BSK19 


17.6312 


16.6696 


12.7398 



as conservative from this point of view. Note also that 

(2) 

even using the value a 2 = does not degrade by more 
than 5% the measurability of Gfi 2 . This is due to the 
presence of a rather large contribution (6.99) in the co- 
efficient of x 2 in Eq. (35) giving the PN correction factor 
to the tidal phase. 



V. CONCLUSIONS 

To conclude, let us summarize the main steps of our 
analysis as well as our main results. 

• For describing the non-tidal, point-mass contribu- 
tion to the Fourier domain phasing, <^°(/), we used 
a 2PN TaylorF2 approximant, which suffices for 
our purposes given the fact that the measurability 
of the nontidal parameters (M, v) is mainly con- 
tained in a rather early part of the inspiral signal, 
say for frequencies / < 50 Hz (see Fig. 3). 

• By contrast, as the tidal contribution to the Fourier 
domain phasing & T (f) is dominated by the late- 
inspiral part of the waveform, wc had to work 
harder and to use the tidal-EOB formalism as a 
way to define a controlled, analytical description 
of the phasing of tidally interacting BNS systems 
up to merger. The controlled analytical description 
we use is a mixture of stationary phase approxima- 
tion and of suitably accurate post-Newtonian ex- 
pansions. We have checked that our approximation 
to the tidal EOB dephasing is accurate to better 
than 0.3 rad up to merger. 

• The final total phasing <&(f) = &°{f) + W T (f) de- 
pends on a number of parameters. We first used 
observational data of known binary pulsars to set 
an apriori upper bound on the magnitude of the 
spin-dependent effects of <?°(/). These effects are 
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proportional to a spin-orbit parameter, (3 and a 
spin-spin one, a. We found that observational data 
suggest that \0\ < 0.2 and \a\ < 10~ 4 . We showed 
that enforcing these upper bounds as Bayesian pri- 
ors in a data analysis is essentially equivalent to 
neglecting from the start /? and a . 

• We have computed G/j.2 for a large sample of EOSs 
(including three recently defined EOSs: SK19, 
BSK20 and BSK21 [43]). 

• The previous considerations allowed us to perform 
a data analysis based on a Fisher matrix formalism 
containing five parameters {t c , </> c , M, v, Gfi2}, and 
taking into account an EOB-controlled analytical 
GW signal going up to merger. 

• The main result of our analysis is that the tidal 
polarizability coefficient G\i2 can be measured at 
the 95% confidence level by the advanced LIGO- 
Virgo detector network using GW signals with rea- 
sonable SNR (p = 16). This mcasurability result 
holds for all the EOSs in the sample we have con- 
sidered, under the only restriction that their maxi- 
mum mass be larger than the recently observed NS 
mass 1.97M© [25]. This measurability property is 
true for BNS mass ratios at least between 0.7 and 
1. 

• We proposed a promising new way of extract- 
ing EOS-dependent information from the coherent 
analysis of a collection of GW observations of sep- 
arate BNS merger events. 

• The latter method is based on parametrizing the 
unknown EOS-dependent function G/j,2(Ma) by a 
(local) linear fit, Eq. (76), depending on only two 
parameters: (oeoSj^eos)- These two parameters 
essentially contain all the information about the 
EOS of neutron star matter that can be extracted 
from GW observations. It would be interesting to 
study the map between (oeoSi^eos) and various 
parametrizations of the EOS introduced in the lit- 
erature [50]. 

• In this paper we have focused on BNS systems, 
but our formalism can be used for discussing mea- 
surability of tidal parameters in mixed black hole- 
neutron star (BHNS) binary systems. However, 
while we have seen that for an equal-mass BNS 
systems the tidal dcphasing at contact was of or- 
der — 10 rad (as analytically given by the approxi- 
mate formula Eq. (32)), for realistic mixed BHNS 
systems this dephasing turns out to be smaller by 
approximately two orders of magnitude. More pre- 
cisely, when considering BHNS systems with mass 
ratio q = Mb /Ma > 1 the tidal dcphasing at con- 
tact (formally defined by Eq. (36)) can be written 
as the product of the equal-mass result (32) by a 



factor Fa, 
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where the supplementary factor Fa is given by 
4 l + 12q 1 



Fa 



13qy/T+q, 1 + qCAC -if^ 



(84) 



with Cb = 1/2 denoting the compactness of the 
black-hole and Ca ~ 1/6 that of the NS. For in- 
stance, for a canonical BHNS system with Ma = 
1AM Q , C a = 1/6 and M B = 1OM , so that 
q ss 7.14, this formula predicts that the BHNS tidal 
dcphasing at contact will be smaller than the BNS 
one by a factor 1/16, so that the BNS, 10-radians 
tidal signal at contact becomes reduced to a level 
~ —0.6 rad. This reduction by more than one order 
of magnitude of the inspiral tidal dephasing signal 
indicates that (even if the coherent analysis sug- 
gested above allows one to work with large effective 
SNRs) it will probably not be possible for LIGO- 
Virgo observations to extract a useful measurement 
of G/j.2 from the inspiral signal. Note that the fac- 
tor Fa is a strongly decreasing function of q. For 
instance, for the extreme mass ratio q = 3 (corre- 
sponding to a 4.2M Q + IAMq system), the factor 
Fa is equal to 0.33, leading to a maximum dephas- 
ing of order —3.3 rad. In addition, we should re- 
member that the measurability of G/12 effectively 
uses only a rather small fraction of the total tidal 
dephasing signal. As discussed in Sec IV D, this 
small fraction is due to the fact that only a small 
part of tf fT (/) is "orthogonal" to the signals associ- 
ated to the nontidal parameters A a - [For instance in 
the equal-mass BNS systems discussed above only 
about 0.1 rad of the maximum 10 rad were use- 
ful in determining the measurability of G\ii, see 
Eq. (75)]. 

• Though our work has been using several simplifica- 
tions (analytical approximations to both the non- 
tidal and tidal part of the EOB Fourier domain 
phasing, Fisher matrix analysis), the various checks 
we have done make us confident that our main 
conclusions are robust. We leave to future work 
a fuller study using better approximations (direct 
Fourier transform of the full, numerically-generated 
EOB waveform, Monte-Carlo estimate of statistical 
errors,. . . ). 

• Concerning our (binary-pulsar-data based) as- 
sumptions about the relative smallness of the spin 
parameters (j3,a), they will be verifiable (or falsi- 
fiablc) once BNS inspiral signals with sufficiently 
high SNRs become available: both through consis- 
tency checks of various measurements of tidal pa- 
rameters, and through direct (Bayesian) analysis of 
the preferred a posteriori values of (/3, a). 
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Appendix A: Description of tidal effects in the 
effective one body model 

This Appendix is devoted to collect all the technical 
elements that are needed to include the description of 
tidal effects in the EOB formalism. The aim of presenting 
this (partly well known) material here, in a self-contained 
form, is twofold: first, it serves to give the reader the 
necessary background to understand how the EOB curve 
of Fig. 1 was generated; second, it is a compact reminder 
of useful formulas in order to help the interested reader 
to do his own EOB implementation. 



1. Dynamics and waveforms 

The EOB formalism [12, 13, 15] replaces the PN- 
cxpanded two-body interaction Lagrangian (or Hamilto- 
nian) by a resummed Hamiltonian, of a specific form, 
which depends only on the relative position and momen- 
tum of the binary system. For a nonspinning BBH sys- 
tem, it has been shown that its dynamics, up to the 3PN 
level, can be described by the following EOB Hamiltonian 
(in polar coordinates, within the plane of the motion): 



H EO B(r,Pr„P v ) = Mc 2 ^\ + 2v(H cil - 1) 
where 



(Al) 



\A+v^ 



square-root" structure Heob ~ y 1 + \JP 2 + ' ' " an d (ii) 
the "condensation" of most of the nonlinear rclativistic 
gravitational interactions in one function of the (EOB) 
radial variable: the basic "radial potential" A(r). The 
structure of the function A(r) is rather simple at 3PN, 
being given by 



(A3) 



1/r 



A^ N (r) = l-2u + 2vu A + a i vu i 

where a 4 = 94/3 - (41/32)tt 2 , and u 
GM/(c 2 tab)- R was recently found that an excellent de- 
scription of the dynamics of BBH systems is obtained [19] 
by: (i) augmenting the presently computed terms in the 
PN expansion (A3) by additional 4PN and 5PN terms; 
(i) Padc-resumming the corresponding 5PN "Taylor" ex- 
pansion of the A function. In other words, the BBH (or 
"point mass") dynamics is well described by a function 
of the form 



a$vu 5 + a^vu 6 ] 



A°(r) = Pi [l - 2u + 2vu z + a 4 vu 4 

(A4) 
where P^ denotes an (n,m) Pade approximant. It was 
found in Ref. [19] (and then substantially confirmed by 
Rcf. [21]) that a good agreement between EOB and 
numerical-relativity BBH waveforms is obtained in an 
extended "banana- like" region in the (05, a 6 ) plane ap- 
proximately spanning the interval between the points 
(a 5 ,a 6 ) = (0,-20) and {a 5 ,a 6 ) = (-36, +520). In this 
work we will select the values 05 = —6.37, a^ = +50, 
which lie within this region (the use of other values within 
the "good BBH fit" region would have no measurable in- 
fluence on the dynamics in the presence of dynamical 
tidal effects). 

The proposal of Ref. [18] for including dynamical tidal 
effects in the conservative part of the dynamics consists 
in simply using Eqs. (A1)-(A2) with the following tidally- 
augmented radial potential 



A(u)=A {u) + A tldal (u). 



(A5) 



Here A°(u) is the point- mass potential defined in 
Eq. (A4), while A tldal (u) is a supplementary "tidal con- 
tribution" of the form 



fl"«ff = W«a 



! + % + *3^f 



Pr, 



(A2) 



Here M = Ma + Mb is the total mass, v 
Ma Mb /{Ma + Mb) 2 is the symmetric mass ratio, and 
23 = 2^(4— 3z/). In addition, we are using rescalcd dimen- 
sionless (effective) variables, namely r = r abc 2 I GM and 
p v = P v c/ (GMaMb), and p rr is canonically conjugated 
to a "tortoise" modification of r [17]. 

A remarkable feature of the EOB formalism is that the 
complicated, original 3PN Hamiltonian (which contains 
many corrections to the basic Newtonian Hamiltonian 
i 2 — 1/?') can be replaced by the simple structure (Al)- 
(A2), whose two crucial ingredients arc: (i) a "double 
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tidal 
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Kiu M+2 A\ l ^{u) , 



(A6) 



= where the terms nju 2i+2 represent the leading-order 
(LO) tidal interaction, i.e., the Newtonian order tidal 
interaction. 

The additional factor A^ dai (u) in Eq. (A6) represents 
the effect of distance dependent, higher-order relativistic 
contributions to the dynamical tidal interactions: 1PN 
(first order in u, or next-to-lcading order), 2PN (of order 
u 2 , or next-to- next-to-lcading order) etc. Here we will 
consider it written in the Taylor-expanded form 



A? i *\u) = l + a? ) u + a¥ ) u 2 



(A7) 
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where ai' are functions of Ma, Ca, and kf for a general 
binary and are written as (see Eq. (37) of [18]) 



-M 



A A(l) . B B{1) 



(A8) 



AW 



where the a n y ' is the coefficient of the nPN fractional 
correction to the tidal interaction potential of body A. 
(see Sec. IIIC of [18]). The dimensionless coefficients 

a n is a function of the dimensionless ratio Xa = 
Ma/M. The analytical expression of the (£ = 2) coeffi- 
cient a 1 has been reported in [18] (and then confirmed 
in [28]) and reads 



A(2) 



= 2 XA ' 



(A9) 



(21 

which in the equal-mass case, Xa = 1/2, yields a\ = 
1.25. Reference [24] has succeeded in computing the first 
post Newtonian octupolar {£ = 3) coefficient a 1 and 
the second post- Newtonian quadrupolar (£ = 2) and oc- 
tupolar (£ = 3) coefficients a 2 . We recall here only the 
most relevant, 2PN quadrupolar one, that reads 



.,^(2) 
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(A10) 



In the equal-mass case, Xa = 1/2, the values of these 



coefficients arc a 



A(2) 



-12) 
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5/4 = 1.25 and a. 



A(2) 



a 2 ' - 85/14 rj 6.071429. In the main text we have 
considered only the tidal (1PN and 2PN) quadrupolar 
contributions, i.e. we have consider only the £ = 2 value 
in Eqs. (4) and (9). In Sec. A 2 below we will investigate 
the (small) effect of the higher-^ tidal corrections for an 
equal-mass BNS. To do so, we adopt the simplifying as- 
sumption that the higher-multipolar tidal-amplification 
factors ^ ldal (u), for £ > 2, are taken to coincide with 
the £ = 2 one. This means that the EOB model that 
we will use here contains only a 1 and a 2 higher 
order tidal parameters that are taken to replace the var- 
ious ah, , with £ = {2,3,4,...}, entering Eq. (A7), i.e. 
^tidai = ^tidai f pr £ > 2. Since the main effect is due 
to the leading-order, Newtonian prefactor kJu 21+2 , this 
simplifying choice does not change in a relevant manner 
the conclusions of Sec. A 2. 

Now that we have reminded the important elements 
needed to build the tidally extended EOB Hamiltonian, 
let us move to discuss how the point-mass EOB waveform 
h® m is augmented by tidal contributions. Similarly to 
the additive tidal modification (A5) to the A potential, 
we will consider an additive modification of the wave- 
form [23], having the structure 



him — h im + h 



tidal 



(All) 



The point-mass contribution is explicitly given by [16] 



ci+ 



W^S^h^p\ m h^ C , (A12) 



where e := ir(£ + m) = 0, 1 is the parity of the considered 
multipolc, where the coefficient 



c i+e (»)=x^ + (-y+*x{+*-\ 



(A13) 



has been separated off the Newtonian waveform h[' — 
cg +t (y)h^ m , and where the other factors respectively 

represent: a source factor S^\ a tail factor h^, a re- 
summed modulus correction p\ m , and a next-to-quasi- 
circular correction, /i.J/ . The latter correction con- 
tains two next-to-quasi-circular parameters (01,02) as 
in Ref. [19]. [Since we will be dealing with equal-mass 
binaries, we fix a\ = —0.0439 and a 2 = 1.3077, ac- 
cording to the EOB/NR comparison (for a BBH equal- 
mass system) of Ref. [19]]. The tail factor introduced 
here is given by h\^ — Ti m e lStm , according to the no- 
tation of Ref. [16]. Using the recent computation [30] 
of the lPN-accurate Blanchet-Damour mass quadrupole 
moment [51] of a tidally interacting binary system (to- 
gether with the Newtonian-accurate spin quadrupole and 
mass octupole) and transforming their symmetric-trace- 
free tensorial results into our £m-multipolar form, we 
have computed the corresponding lPN-accurate value of 
/i^ dal , as well as the OPN-accurate values of h 2 f al , ft^ dal , 
and hqi- In addition, using the general analysis of tail 
effects in Refs. [52, 53] and the resummation of tails in- 
troduced in Refs. [16, 54, 55], we were able to further 
improve the accuracy of these waveforms by incorporat- 
ing (in a resummed manner) the effect of tails (to all 
orders in M). From a PN point of view, this means, in 
particular, that the tidal contribution we use to the to- 
tal metric waveform is 1.5PN accurate. From the results 
of Ref. [30], 1PN source moment, one has that the only 
nonvanishing multipolar tidal corrections at 1PN frac- 



(m = 1,2) and /^ ldal ( 



1,3). 



tional level are /4m al 
These 1 multipolar components of the waveform can be 
obtained simply by computing the corresponding num- 
ber of time derivatives (in the circular approximation) 
of the multipole moments given by [30] and then pro- 
jecting them along symmetric-trace-free tensor spherical 
harmonics. For consistency with the fact that tidal effects 
are included at 2PN fractional accuracy in the Hamilto- 
nian, we similarly write the (multipolar) waveform such 
to formally include (yet analytically uknown) 2PN tidal 
corrections. The £ = m = 2 tidal part of the waveform is 
written as 



/i" daI = h 



r(N,0) 
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v™ll + ft\XB)vl + P?(X B y i2 \ 



(A14) 



where the f3 22 (X) and (3 22 (X) functions parametrize re- 
spectively 1PN and 2PN fractional tidal corrections. The 
1PN tidal coefficient, /3 22 (X), can be computed analyti- 
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cally from the results of Ref. [30] and it reads 



Pf(X) 



-202 + 560A - 340A 2 + 45A E 
42(3 - 2X) 



(A15) 



The 2PN tidal coefficient /3f 2 (A) is currently unknown 
analytically. Note that, following the original suggestion 
of Ref. [56], in Eq. (A14) we have replaced the PN order- 
ing parameter x = (A/ J?) 2 / 3 by the EOB velocity variable 
vq = ra£2, where ra = rip 1 / 3 and ip is computed using 
the 3PN-accurate EOB Hamiltonian following the defi- 
nition (originally at 2PN accuracy) of Ref. [56] (see also 
Ref. [17]). The other waveform multipoles that present 
tidal corrections (up to the 2PN formal level) are 

/i 2 i = h' 2 [ ' '< c 2 +i(y)h 2 i 



--6X B )vjf[l + ^ 1 (X B )vj i }Ki 



"--&X A )v 1 £[l + Pl 1 {X A )vl]4 ( 



h 3m = ^m <i c 3+0 {v)h 3m 



+ 6X A v 1 a °(l + {3 3m (X B )v 2 n )K 



(A16) 



6A B ^°(l+^ 3m (X A K) K 



(A17) 



where /3f m (A) formally indicate the (currently unknown) 
corresponding 2PN corrections. The tidally-corrccted ra- 
diation reaction J- v is then computed by using this cor- 
rected waveform in the definition of T v given in [16, 19]. 
Note that the point-mass partial multipolar amplitudes 
pe m that we use in the construction of the analytic radi- 
ation reaction in this paper are augmented with respect 
to those discussed in Ref. [16, 19] by the new (5PN ac- 
curate) v = terms recently computed in [57]. By con- 
trast, for simplicity and for consistency with Ref. [19], 
we adopt the same prescription of that reference, Eq. (4) 
there, to compute the i — m — 2 point-mass wave- 
form, which relies on a different resummation of the 
residual amplitude correction with respect to Eq. (A12) 
above. More precisely, the residual modulus correction 
of Eq. (A12), (P22) 2 , is replaced by the Pade-resummed 
function fH(x;v) = Pitf^^ix; v)], where / 2 1 2 aylor is 
computed in Ref. [16] at 3 +2 PN accuracy. In addition, 
the residual phase correction 822 that is used here is com- 
puted at the accuracy given in [16], Eq. (20), without the 
further (y — 0) term obtained by [57]. 

In summary, the EOB tidal model that we use here 
is formally complete up to the 2.5PN level, though an- 
alytically complete at the 1.5PN level only, because of 
the (current) lack of analytical information on the coeffi- 
cients {Pf{X), /3 21 (A), /? 31 (A), /?f (X)}. Despite this, 
one should keep in mind that in the most relevant equal- 
mass case, only the /3| 2 coefficient is relevant. On top of 
this, let us remind that in Sec. II C we argued that in the 
(Fourier domain) tidal phasing the contribution due to 
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FIG. 5. Difference between the full SPA tidal phase ^ 4 
obtained from a EOB waveform computed with an ^4 tldal in- 
cluding &2, kz and ki, and the tidal phase computed with ki 
only. The figure refers to a 7 = 2, rest-mass polytrope BNS 
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model with C = 0.16 and a.__ 
indicates the contact frequency. 



6. The vertical dashed line 



/3| 2 is very subdominant with respect to the others and 
thus it can be safely neglected in first approximation. 

In the main text we have used the EOB approach 
to compute a SPA-dcfined, Fourier domain EOB tidal 
phasing, !?eoBsp A ( w )' where w is the quadrupolar GW 
frequency so as to control the accuracy of the Fourier- 
domain, PN-expanded tidal phase given by Eq. (31) for 
an equal-mass binary. The analytical procedure to ob- 
tain ^eoBspa^) (t na t we shall henceforth simply de- 
note as }& T {<jj)) is described in Sec. II B (see in particular 
Eq. (20)), while the comparison with the PN-expanded 
tidal phasing is discussed in Sec. (II C), in particular 
Figs. 1 and 2 there. The phase ^ T (w) is computed by 
integrating numerically Eq. (19) starting from the fre- 
quency ujq that marks the beginning of the inspiral wave- 
form obtained when solving the EOB equations of motion 
numerically. This integration is done using the 2.5PN re- 
sult for $^ 5 pN and dW^ 5PN / ' dui as initial boundary con- 
ditions, and thus loq needs to be chosen sufficiently small 
(i.e., the iniitla separation is sufficiently large) so to have 
Q2.5PN py qEOB^ j n p rac tice, for all compactnesses con- 
sidered in Figs. 1 and 2, the initial relative separation 
that we use is r$ — 32, which corresponds to a quadrupo- 
lar GW frequency uj ~ 0.111. 



2. Dependence of the phasing on dynamical tidal 
effects with £ > 2 and its linearity with K2 

The use of the tidal phasing ^spn- Eq. (31) in the 
main text neglects by construction two physical effects 
that are incorporated in the EOB description, notably: 
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\Pk computed when retaining only &2- One sees that the 
phase difference at the contact frequency (indicated by 
the dashed vertical line) is of the order of —0.2 rad. 

Finally, to study to what extent the dependence of &£ 
is linear in &2 (and therefore K2), we consider three bina- 
ries with compactnesses C = (0.14,0.16,0.18). For each 
binary we compute the ratio W^ Jk^, that is displayed in 
Fig. 6. This hgure indicates that &£ [ s linear in k^ to 
good approximation. From the difference (—0.0015 rad) 
between the C = 0.18 and C = 0.14 curve at Mtu = 0.1048 
(contact of C = 0.14 binary) one can estimate that tak- 
ing into account the kJ would further decrease \]/£ by no 
more than —0.3 rad. Again, neglecting this effects means 
that our estimates will be on the conservative side. 

Appendix B: The tidal phase for a general binary at 
2.5PN accuracy 



FIG. 6. Approximate linearity of $ec>b with respect to k 2 , 
with aj =6. Three equal-mass (polytropic) BNS systems 
of compactnesses C — {0.14, 0.16, 0.18} are compared. The 
vertical dashed lines indicate the corresponding contact fre- 
quencies. 



(i) dynamical tidal terms with £ > 2 (e.g., £ = 3 and 
£ = 4) that enter in the definition of A tldal (u); (ii) non- 
linear effects in /c^, that are present due to the resummed 
nature of the EOB formalism. In this Section we analyze 
their influence of !? T (u;) and argue that they collectively 
contribute to the phasing by an amount of the order of 
—0.2 rad up to contact. On top of this contribution being 
small, the fact that it is negative means that any measur- 
ability analysis neglecting it is on the conservative side. 
Focusing on a C — 0.16, rest-mass, 7 = 2, polytrope bi- 
nary, the effect of the I = 3 and I = 4 tidal corrections 
to the EOB potential is illustrated in Fig. 5. The fig- 
ure shows the difference between the full \PT obtained 

^2,3,4 

from a EOB waveform computed retaining &2, &3 and 
&4 in the Hamiltonian ( through A tldal ), and the phase 



In this Appendix we collect PN-cxpanded expressions 
for Q^ and for ^ T (w) for a general binary. Such a result 
is here expressed as a function of the PN ordering param- 
eter x, of Xa = Ma/M and Xb = 1 — Xa and of the di- 
mensionless tidal parameter k^ defined in Eq. (6). Note 
that this result is general, in the sense that it holds for a 
neutron star binary of any mass ratio, or a mixed, black- 
hole neutron star binary. In this latter case the tidal 
parameter of one of the two objects is put to zero [6]. 
The equal-mass analytic expressions used in the main 
text are obtained as a particular case of the equations 
listed below, i.e., Xa = Xb = 1/2. 

In detail, the 2.5PN accurate tidal part of the Q u (x) 
function, Q^(x), has the structure 



)l = Qt{x)i4 



\{x)4 



(Bl) 



where the Q^(x) contribution is written as the following 
PN series 

Qt = 9Ncwt^ 5/2 (l + fix + q£x 3/2 + q£x 2 + q£x 5/2 ) , 

(B2) 
where the coefficients read 
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X B 



5 12 - 11X A 



3179 - 919X A - 2286X5 + 260X5 



24z/ 



X A 

3 

.•1 



48(12 -IIXa) 



-4-7T, 



(B3) 



(B4) 
(B5) 



1 



VLX A - 12 



67702048X5 - 223216640X;, + 337457524X5 - 141992280X5 + 96008669X A - 143740242 



338688 



( l A X -\x 2 + ^- 1 -) 



(1-2X j4 )(1-X j4 ) 2 



Pl l {X A ) + ^l\Xs 



it 



1344'" L v "~' ' 448 
7T (27719 - 22127X A + 7022X^ - 10232X 



96(12 -11X A ) 



(B6) 

(B7) 



r 



and we left explicit the dependence on the (yet unknown) 
parameters entering the tidal contribution to the wave- 
form at fractional 2PN order, {/3f 2 , /? 21 , f3f\ (if). The 
tidal part of the phase is then written in the following 
form 



<£ 



2.5PN 



*H?)i4 + *b(x)k 



LI 
2 • 



(B8) 



where the 'F^(x) contribution is given by the following 
PN series 



iZ/ T - TV 4 



x 5/2 (l + pfx + p2X 3/2 + pfx 2 + p£x 5/2 ) , 



whose coefficients read 
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PNowt 



I6v 



12 



Xa 
X b 



3 12 - 11X A 
V&v 1-X A '■ 



_ 5 (3179 - 919X^4 - 2286X A + 260X A ) 



Pi 






672(12 -UX A ) 
1 f 39927845 480043345 



12 - \\X A 

5 



508032 



9144576 



X A 



9860575 
127008 



X 



421821905 



2286144 



X 



2 -X A \ pf + A (3 - 13X A + 18X 2 - 8Xi) /3 21 (X A ) 



1 



+ (1 - X A f (1 - 2X A ) 



5588^™ + ^™ 



ft 1 



tt(27719 - 22127Xa + 7022X^ - 10232X^) 
672(12-11X^4) 



4359700 
35721 



X 4 



10578445 
285768 



(B10) 



(Bll) 
(B12) 



X 



(B13) 



r 



The 1PN coefficient pf above was obtained in 
Eqs. (3.8a)-(3.8b) of [30]. Then the relation 



,GtijX B 

' M 5 X A ' 



(B14) 



(and similarly with i f) B) is used so to express the tidal 
phase in terms of the tidal polarizability coefficients GpJ^ 



and Gfif ■ By defining 



GU 2 = — 
^ 26 



1 



12 f » 



1 



12^) 



(B15) 
a straightforward calculation finally yields the tidal part 
of the phase in the form 



^2.5PN 



8vM 5 



2.5PN 



(B16) 
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that is the one used (with 



3? 



ft 1 = PI 3 = 0? = 0) 



in Eq. (60) to perform the Fisher-matrix-based estimate 
of the measurability of G/I2 • 
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